New approach to the treatment of separatrix chaos 
and its application to the global chaos onset between adjacent separatrices 
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We have developed the general method for the description of separatrix chaos, basing on the 
analysis of the separatrix map dynamics. Matching it with the resonant Hamiltonian analysis, we 
show that, for a given amplitude of perturbation, the maximum width of the chaotic layer in energy 
may be much larger than it was assumed before. We apply the above theory to explain the drastic 
facilitation of global chaos onset in time-periodically perturbed Hamiltonian systems possessing 
two or more separatrices, previously discovered (PRL 90, 174101 (2003)). The theory well agrees 
with simulations. We also discuss generalizations and applications. Examples of applications of 
the facilitation include: the increase of the DC conductivity in spatially periodic structures, the 
reduction of activation barriers for noise-induced transitions and the related acceleration of spatial 
diffusion, the facilitation of the stochastic web formation in a wave-driven or kicked oscillator. 



PACS numbers: 05.45.-a, 05. 45. Ac, 05.45.Pq 
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I. INTRODUCTION 

A weak perturbation of a Hamiltonian system causes 
the onset of chaotic layers around separatrices of the un- 
perturbed system and/or separatrices surrounding non- 
linear resonances generated by the perturbation 0, 0, 0, 
0, 0]. The system may be transported along the layer 
in a random-like fashion and this chaotic transport pkws 
an important role in many physical phenomena [H, 0, 0] . 
If the perturbation is sufficiently weak, then the layers 
are thin and the chaos is called local [l|, 0, 0, 0] . As the 
perturbation magnitude increases, the width of the layer 
grows and the layers corresponding to adjacent separatri- 
ces reconnect at some, typically non-small, critical value 
of the perturbation. This conventionally marks the onset 
of global chaos 0, 0, 0, 0] he. chaos in a large region of 
the phase space, with chaotic transport throughout the 
whole relevant energy range. 

The reconnection of the layers around separatrices of 
the resonances often correlates with the overlap in en- 
ergy between neighbouring resonances calculated inde- 
pendently in the resonant approximation. The latter 
constitutes the heuristic Chirikov resonance-overlap cri- 
terion 0, 0, 0, 0] . Rut the Chirikov criterion may fail if 
the system is of the zero- dispersion (ZD) type [6| i.e. if 
the frequency of eigenoscillations possesses a local max- 
imum or minimum as a function of its energy (cf. also 
studies of related maps 0, 0] which are called nontwist, 
twistless or nonmonotonic twist maps). In such systems, 
there are typically two resonances of one and the same 
order [9( , and their overlap in energy does not result in 
the onset of global chaos [a 0, HI • Even their overlap in 
phase space [10| results typically only in the reconnection 




*Also at Physics Department, Lancaster University, UK 



Slow an^le 



FIG. 1: Typical evolution of thin chaotic layers in the plane of 
slow variables of a zero-dispersion system (the perturbation 
magnitude grows from the top to the bottom). 



of the thin chaotic layers associated with the resonances. 
As the amplitude of the time-periodic perturbation grows 
further, the layers may separate again |,0,|. An exam- 
ple of the evolution of resonances in the plane of energy 
and slow angle is given in Fig. 1 (the typical evolution 
of a real Poincare section is shown e.g. in [111]). 

As it is known Q , any Hamiltonian system with two or 
more separatrices belongs to the ZD type: the eigenfre- 
quency as a function of energy possesses a local maximum 
between each pair of adjacent separatrices. For the pur- 
pose of global chaos onset, our letter [l^] has addressed 
the possibility to combine the overlap of resonances with 
each other (typical of ZD systems) and their overlap with 
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the chaotic layers associated with the separatrices. Via 
numerical simulations, [l2| demonstrated that this is pos- 
sible, leading to a scenario for global chaos onset which 
requires much smaller perturbation amplitudes than in 
the conventional case. The letter [12J suggested also a 
heuristic theory for this effect (more details were pre- 
sented in [13|). 

The present work develops the method for the quan- 
titative description of chaotic layers in phase space, for 
the resonance frequency range. We uncover the physical 
mechanism of their overlap with the resonances, and on 
this basis develop a detailed self-contained theory of the 
facilitated onset of global chaos. We also discuss general- 
izations and applications. The method for the description 
of the chaotic layers is of general importance: in conven- 
tional systems with a single-separatrix layer, it predicts 
a much larger maximum width in energy than what was 
assumed before [H, 1H|. 

The paper is organized as follows. Sec. II introduces 
some relevant model example and presents the major re- 
sults of the simulations: studying numerically the fre- 
quency dependence of the minimal amplitude of the AC 
drive for which global chaos occurs, h gc {ojf), we show 
that h gc (ujf) possesses deep spikes at certain frequencies. 
Sec. Ill gives the self-consistent asymptotic theory for 
the minima of the spikes, after assessing the boundaries 
of the relevant chaotic layers. Sec. IV gives the theory 
for the spikes wings. Discussion of a few generalizations 
and applications is carried out in Sec. V. Conclusions are 
drawn in Sec. VI. The Appendix describes in details the 
new method for the analysis of separatrix chaos. 



II. MODEL AND MAJOR RESULTS OF 
SIMULATIONS 

As an example of a one-dimensional Hamiltonian sys- 
tem possessing two or more separatrices, we use a spa- 
tially periodic potential system with two different-height 
barriers per period (Fig. 2(a)): 
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This model may relate e.g. to a pendulum spinning 
about its vertical axis [l4| or to a classical 2D electron 
gas in a magnetic field spatially periodic in one of the 
in-plane dimensions [HI, TiH ] . The interest to the lat- 
ter system arose in the 90th due to technological ad- 
vances allowing to manufacture magnetic superlattices 
of high-quality [13, EH leading to a variety of interest- 
ing behaviours of the charge carriers in semiconductors 
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Figs. 2(b) and 2(c) show respectively the separatrices 
of the Hamiltonian system (1) in the p— q plane and the 
dependence of the frequency lj of its oscillation, often 
called eigenfrequency, on its energy E = Hq (p, q) . The 
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FIG. 2: The potential U(q), the separatrices in the phase 
space and the eigenfrequency cu(E) for the unperturbed sys- 
tem (1) with <E> = 0.2, in (a), (b) and (c) respectively. 



separatrices correspond to energies equal to the value of 
the potential barrier tops E^ = (1 — $) 2 /2 and E^ 2 ' = 
(1 + $) 2 /2 (Fig. 2(a)). The function lo{E) is close to 
the extreme eigenfrequency u/ m = uj{E m ) for most of the 

range [E^ 1 , E^ 2 ] while sharply decreasing to zero as E 

approaches either or E^ . 

Add now a time-periodic perturbation: as an exam ple , 
we use an AC drive, which corresponds to a dipole [3.l2lj 
perturbation of the Hamiltonian: 

q = dH/dp, p = -dH/dq, (2) 
H(p, q) = H (p, q) - hqcos(ujft). 

The conventional scenario of global chaos onset be- 
tween the separatrices of the system (2)-(l) is illustrated 
by Fig. 3. The figure presents the evolution of the strobo- 
scopic Poincare section as h grows while loj is fixed at an 
arbitrarily chosen value away from oj m and its harmon- 
ics. At small h, there are two thin chaotic layers around 
the inner and outer separatrices of the unperturbed sys- 
tem. Unbounded chaotic transport takes place only in 
the outer chaotic layer i.e. in a narrow energy range. As 
h grows, so do also the layers. At some critical value 
h gc = h gc (ujf), the layers merge. This may be consid- 
ered as the onset of global chaos: the whole range of 
energies between the barrier levels is involved, with un- 
bounded chaotic transport. The states {1^} = {p = 
0,q = tt/2 + 2ttZ} and {0«} = {p = 0, q = -tt/2 + 2ttZ} 
(where I is any integer) in the Poincare section are associ- 
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FIG. 3: The evolution of the stroboscopic (at t = n2n/uf 
with n = 0, 1, 2, ...) Poincare section of the system (2)-(l) with 
$ = 0.2 as h grows while uJf = 0.3. The number of points in 
each trajectory is 2000. In (a) and (b), three characteristic 
trajectories are shown: the inner trajectory starts from the 
state } = {p — 0, q — 7r/2} and is chaotic but bounded in 
space; the outer trajectory starts from {O^} = {p = 0,q = 
— vr/2} and is chaotic and unbounded in coordinate; the third 
trajectory is an example of a regular trajectory separating the 
two chaotic ones. In (c), the chaotic trajectories mix. 

ated respectively with the inner and outer saddles of the 
unperturbed system, and necessarily belong to the inner 
and outer chaotic layers, respectively. Thus, the neces- 
sary and sufficient condition for global chaos onset may 
be formulated as the possibility for the system placed ini- 
tially in the state {1^} to pass beyond the neighbouring 
of the "outer" states, {O (0) } or {0 (1) }, i.e. the coordi- 
nate q becomes < — ir/2 or > 3tt/2 at sufficiently large 
times t 3> 2n/ujf. 

A diagram in the h — ujf plane, based on the above 
criterion, is shown in Fig. 4. The lower boundary of 
the shaded area represents the function h gc (uif). It has 
deep spikes i.e. cusp-like local minima. The most pro- 
nounced spikes are situated at frequencies uif — to/ that 
are slightly less than the odd multiples of uj m , 

« w m (2j - 1), j — 1, 2, ... (3) 

The deepest minimum occurs at uffl ~ u) m : the value of 
h gc in the minimum, ha = h gc (u>i ), is approximately 
40 times smaller than the value in the neighbouring pro- 
nounced local maximum of h gc (ujf) at w/ ~ 1. As n 
increases, the nth minimum becomes less deep. The func- 
tion h gc ((jjf) is very sensitive to w/ in the vicinity of the 

minima: for example, a shift of ojf down from u>i ~ 0.4 
by only 1% causes an increase of h gc by ~ 30%. 

The origin of the spikes becomes more clear looking 
at the evolution of the Poincare section for Wf ~ wi 1 ' as 
h grows (Fig. 5): it drastically differs from the con- 
ventional evolution shown in Fig. 3. For h — 0.001 




FIG. 4: The diagram indicating (shade) the perturbation pa- 
rameters range for which global chaos exists. Integration time 
for each point of the grid is 120007T. 

(Fig. 5(a)), one can see four chaotic trajectories. Two 
of them are associated with the inner and outer separa- 
trices of the unperturbed system, similarly to the conven- 
tional case (cf. Fig. 3). They are marked by green and 
blue respectively. These trajectories fill the correspond- 
ing chaotic layers, which will be referred below as the 
"inner" and "outer" separatrix layers respectively. The 
other two chaotic trajectories marked by red and cyan 
are associated with the two nonlinear resonances of the 
1st order. Examples of non-chaotic trajectories separat- 
ing the chaotic ones are shown in brown. As the pertur- 
bation amplitude h increases, the outer separatrix layer 
sequentially absorbs other chaotic trajectories while large 
stability islands (associated with the resonances) arise in 
the layer. At h — 0.003, it has absorbed the red trajec- 
tory: the resulting chaotic layer is shown in blue in Fig. 
5(b). At h — 0.00475, this chaotic layer has absorbed 
the cyan chaotic trajectory: the resulting chaotic layer is 
shown in blue in Fig. 5(c) 22]. Finally, at h = 0.0055 
the latter blue layer has merged with the inner separatrix 
layer [23| (see Fig. 5(d)), i.e. the onset of global chaos as 
defined above has occurred. 

Even prior to the theoretical analysis, one can draw a 
few conclusions from the evolution. Namely, if ujf is close 
to the minimum of the spike of h gc (ujf), then 

1) the onset of global chaos occurs due to the com- 
bination of the overlap of chaotic layers associated 
with nonlinear resonances with each other and the 
overlap of the latter layers with the inner and outer 
separatrix layers; 

2) the width of the nonlinear resonances are large al- 
ready at quite small amplitudes of the perturba- 
tion, so that the overlap with the chaotic layers 
around the original separatrices occurs at unusu- 
ally small perturbation amplitudes; 

3) the onset of the overlap of at least one of the non- 
linear resonances with the outer separatrix layer 
occurs at values of h which are a few times smaller 
than those required for the onset of the overlap with 
the inner separatrix layer. 
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FIG. 5: The evolution of the stroboscopic Poincare section of 
the system (2)-(l) with $ = 0.2 as the amplitude of the per- 
turbation h grows while the frequency is fixed at ujf — 0.401. 
The number of points in each trajectory is 2000. The chaotic 
trajectories starting from the states {1^} and {O*- '} are 
drawn in green and blue respectively. The stable stationary 
points of Eq. (14) (the lst-order nonlinear resonances) are 
indicated by the red and cyan crosses. The chaotic layers as- 
sociated with the resonances are indicated in red and cyan 
respectively, unless they merge with those associated with 
the green/blue chaotic trajectories. Examples of regular tra- 
jectories embracing the state {1^} while separating various 
chaotic trajectories are shown in brown. 

The above conclusions are also illustrated by Fig. 6 
which presents the evolution of the phase space of slow 
variables l&llll, action / = 1(E) and slow angle 
ip = ip — lu ft, calculated in resonance approximation for 
the lst-order spike (see Eq. (4) below). 

Similarly, for the spikes of higher order, higher-order 
resonances are relevant. 



III. EXPLICIT ASYMPTOTIC THEORY FOR 
THE MINIMA OF THE SPIKES 

The eigenfrequency lu(E) is close to its local maximum 

Lj m for most of the relevant range [E^\e^] (Fig. 2(c)). 
As shown below, uj(E) approaches a rectangular form in 
the asymptotic limit $ — > 0. Hence, if the perturba- 
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FIG. 6: The evolution of the separatrices (full lines) of the 
lst-order resonances within the resonance approximation (de- 
scribed by (4) with n = 1) in the plane of action I and slow 
angle ip, for the same parameters as in Fig. 5. Horizontal 
levels mark the values of / corresponding to the barriers. 



tion frequency uif is close to Lu m or its odd multiples, 
\trff — (2j — l)w m | <C ui m , then the energy width of non- 
linear resonances becomes comparable to the width of 
the whole range between barriers (i.e. E^ — E^ w 2<f>) 
at a rather small perturbation magnitude h <§C <&. Note 
that $ determines the characteristic magnitude of the 
perturbation required for the conventional overlap of the 
separatrix chaotic layers, when w/ is not close to any odd 

multiple of oj m (Fig. 3 (c)). Thus, if ujf w u)s , the non- 
linear resonances should play a crucial role in the onset 
of global chaos (cf. Fig. 5). 

We note that it is not entirely obvious a priori whether 
it is indeed possible to calculate ha = h gc (uji^) within 
the resonance approximation: in fact, it is essential for 
the separatrices of nonlinear resonances to nearly touch 
the barriers levels, but the resonance approximation is 
obviously invalid in the close vicinity of the barriers; 
furthermore, numerical calculations of resonances show 
that, if cjf w u>s , the perturbation amplitude h at which 
the resonance separatrix touches a given energy level in 
the close vicinity of the barriers is very sensitive to w/, 

apparently making the calculation of hi within the res- 
onance approximation even less feasible. 

Nevertheless, we show below in a self-consistent man- 
ner that, in the asymptotic limit $ — > 0, the relevant 
boundaries of the chaotic layers lie in the range of en- 
ergies E where to(E) w u) m . Therefore, the resonant 
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approximation is valid and it allows to obtain explicit 
asymptotic expressions both for u s and 
the wings of the spikes in the vicinities of u>. 

The asymptotic limit $ — * is the most interesting one 
from a theoretical point of view since this limit leads to 
the strongest facilitation of the global chaos onset and it 
is most accurately described by the self-contained theory. 
Most of the theory presented below assumes this limit 
and concentrates therefore on the results to the lowest 
order in the small parameter. 

On the applications side, the range of moderately small 
$ is more interesting, since the chaos facilitation is still 
pronounced (and still described by the asymptotic the- 
ory) while the area of chaos between the separatrices is 
not too small (comparable with the area inside the inner 
separatrix): cf. Figs. 2, 3 and 5. To increase the accu- 
racy of the theoretical description in this range, we esti- 
mate the next-order corrections and develop an efficient 
numerical procedure allowing for the further corrections. 



tions of motion (2) with Hq (1) can be integrated [ji 
(see also Eq. (60) below), so that we can find u(E): 



w(E) 




(5) 



where 



KM = 



1 - k 2 sin 2 (0) 



(6) 



is the full elliptic integral of the first order |25j . Using its 
asymptotic expression, 



A. Resonant Hamiltonian and related quantities we derive ^(E) in the asymptotic limit $ -> 0: 



Let LJf be close to the nth odd [2J| harmonic of LJ m , n = 
(2j — 1). Over most of the range [eP^E?'], except in 

the close vicinities of E^ and E^\ the nth harmonic of 
eigenoscillation is nearly resonant with the perturbation. 
Due to this, the (slow) dynamics of the action / = 1(E) = 
(27T)" 1 § dqp and the angle if; 0, 1 S 1 S, I, 0, M, E| 
can be shown to be described by the following auxiliary 
Hamiltonian (cf . [j 0, |, 1, 0, g 0, 1 ) : 



u{E) 



In 



(>4 



($-A£)($+AE) 



AE^E- 1 -, 



\AE\ < 



(7) 



$ -> 0. 

The function uj(E) (7) is close to its maximum 



H{I, ip) = / dl (riLj — luj) — nhq n cos(ip) (4) 

Jl(E m ) 

= n(E - E m ) - Uf(I - I(E m )) - nhq n cos(^) , 
I = I(E) = f E E^H (p,q), 

J E, 



ip = nip — LUft, 



ip = it + sgn(p)u(E) 



dq 



qtaln{ E) y/2(E-U(q)) 



2ttL 



2 r /2 

q n = q n (E) = - / dip q(E, ip) cos(n^), 



\nuj — Uf\ <C oj, 



n = 2j-l, j = 1,2,3,.. 



where i? m in is the minimal energy (over all q,p) E = 
Ho(p,q); u> = u(E) = dH^/dl and q m in(E) are, respec- 
tively, the frequency and the minimal coordinate of the 
conservative motion with a given value of energy E; I is 
the number of right turning points in the trajectory [q(r)] 
of the conservative motion with energy E and given ini- 
tial state (qo,po). 

Let us derive the explicit expressions for various quan- 
tities in (4). In the unperturbed case (h = 0), the equa- 



W(E)} 



21n(8/$) 



(8) 



for most of the interbarrier [26[ range of energies [1/2 — 
$, 1/2 + $]; on the other hand, in the close vicin- 
ity of the barriers, where either |ln(l/(l — AE/$>))\ or 
|ln(l/(l + AE/$>))\ become comparable with, or larger 
than, ln(8/<&), uj(E) sharply decreases to zero as \ AE\ — > 
<i>. The range where this takes place is ~ $ 2 , and its ratio 
to the whole interbarrier range, 2$, is ~ $ i.e. it goes to 
zero in the asymptotic limit $ — > 0: in other words, lo(E) 
approaches a rectangular form. As it will be clear from 
the following, it is this almost rectangular form of 
uj(E) which determines many of the characteristic 
features of the global chaos onset in systems with 
two or more separatrices. 

One more quantity which strongly affects (cj s ,h s ) is 
the Fourier harmonic q n = q n (E). The system stays 
most of the time very close to one of the barriers. Con- 
sider the motion within one of the periods of the poten- 
tial U(q), between neighboring upper barriers [q^ b , q^ b ] 
2tt. If the energy E = 1/2 + AE lies in 
then the system will stay 
close to the lower barrier qn, = qty + it for a time [27j 



where q^ b 



a {1) - 



the relevant range [E^\E^ 
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Ti w 21n | 5—) (9) 

during each period of cigcnoscillation, while it will stay 

close to one of the upper barriers q ub ' = qib ± 7r for most 
of the remaining of the eigenoscillation, 



We emphasize that \q n \ determines the "strength" 
of the nonlinear resonances: therefore, apart from the 
nearly rectangular form of lu(E), the non-smallness of 
\q n \ is one more factor giving rise to the strong 
facilitation of the global chaos onset. 

We shall need also the asymptotic expression of the 
action /. Substituting io(E) (7) into the definition of 
1(E) (4) and carrying out the integration, we obtain 



2 In 



1 



$ — AE 



(10) 



Hence, the function q(E, ip) — qib may be approximated 
by the following piecewise even periodic function: 

a(F ,IA n -I* ^ ^4°'^'^M W -5TFfe' 7r ]' 
qy£j,ip) qw — S . Tu Tm , 

q(E, -ip) - qw = q(E, ip) - 3/6, 

q(E,il>±2m,) = q(E,il>), 
i = 1,2,3,... 

Substituting the above approximation for q(E,ip) into 
the definition of g„ (4) , one can obtain: 



taj-i = Q2j-i(E) 



2.7 - 1 



■ sm 



(2j - 1)tt/2 



$ - 

3 = 



>0, 
= 0, 

1,2,3,... 



(11) 



At barrier energies, qij-\ takes the values 



(2.7-1)- 

As i? varies in between the barrier values, Q2j-i varies 
monotonously if j = 1 and non-monotonously otherwise 
(cf. Fig. 11). But in any case, the significant variations 
occur mostly in the close vicinity of the barrier energies 
E^ and E^ while, for most of the range [E9~\E^], 
the argument of the sine in Eq. (11) is close to 7r/4 and 
(fej-i is then almost constant: 



Q2j-1 



(-1)1 



2J-1 1 V2 



2j-l 



1,2,3, (12) 



In 



1 + AE/<f>\ 



1 



«21n(- 



where [. . .] means the integer part. 

In the asymptotic limit $ — > 0, the range of AE 
where the approximate equality (12) for (fcj-i is valid 
approaches the whole range ] — $[. 



AS In 



1(E) = 1(1/2) + 
$ -> 0. 



/ 64e 2 
^$ 2 -(AB) 2 



<f>ln 



(13) 



B. Reconnection of resonance separatrices 

We now turn to the analysis of the phase space of 
the resonance Hamiltonian (4). The evolution of the 
Poincare section (see Fig. 5 and the related analysis 
in Sec. II) suggests that we need to find such separa- 
trix of (4) which undergoes the following evolution as h 
grows: for sufficiently small h, the separatrix does not 
overlap chaotic layers associated with the barriers while, 
for h > h gc (u!f), it does overlap them. The relevance of 
such a condition will be justified further. 

For u)f w nu) m with a given odd n, the equations of 
motion of the system (4) read as follows 



i = -— T = -nhq n sin(r/>), 
dip 

7 dH dq n ~ 

ip = =noj — u>f — n ' l ~^j C0S W)- 



(14) 



Any separatrix necessarily includes one or more unsta- 
ble stationary points. The system (14) may have several 
stationary points per 2-7T interval of ip. Let us first ex- 
clude those points which are irrelevant to the separatrix 
undergoing the evolution described above. 

Given that q n (E\p) = 0, there are two unstable sta- 
tionary points with I corresponding to E = E^ 1 ' and 
ip = ±7r/2. They are irrelevant since, even for an in- 
finitely small h, each of them necessarily lies inside the 
corresponding barrier chaotic layer. 

If E 7^ E^\ then q n ^= 0, so I = only if ip is equal 
either to or to n. Substituting these values into the 

second equation of (14) and putting ip = 0, we obtain 
the equations for the corresponding actions: 



X T (T) = nu — uif =F nhdq n /dl = 0, 



(15) 



where the signs "-" and "+" correspond to ip = and 
ip = it respectively. A typical example of the graphic 
solution of equations (15) for n = 1 is shown in Fig. 7. 
Two of the roots corresponding to ip — ir are very close 
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FIG. 7: A schematic example illustrating the graphic solu- 
tions of Eqs. (15) for n = 1, as intersections of the curve ui(I) 
(thick solid red line) with the curves oj/ ± hdq n (I) / dl (thin 
solid green lines). The solutions corresponding to the lower 
and upper relevant saddles (defined by Eq.(16)) are marked 
by dots and by the labels si and su respectively (we do not 
mark other solutions because they are irrelevant). 



to the barrier values of / (we remind that the relevant 
values of h are small) . These roots arise due to the diver- 
gence of dq/dl as / approaches any of the barrier values. 
The lower/upper root corresponds to a stable/unstable 
point. However, for any n, both these points and the sep- 
aratrix generated by the unstable point necessarily lie in 
the ranges covered by the barrier chaotic layers. There- 
fore, they are also irrelevant [28] . For n > 1, the number 
of the roots of (15) in the vicinity of the barriers may be 
larger (due to the oscillations of the modulus and sign of 
dq n /dl in the vicinity of the barriers) but they all are 
irrelevant for the same reason, at least to leading-order 
terms in the expressions for the spikes minima. 

Consider the stationary points corresponding to the re- 
maining four roots of equations (15). As follows from the 
analysis of equations (14) linearized near the stationary 

points (cf. HUB II all), 

two of them are stable elliptic 
points [29( | . often called nonlinear resonances, while two 
others are unstable hyperbolic points, often called sad- 
dles. These saddles are of main interest in the context 
of our work. They belong to the separatrices separating 
the regions of oscillations around the resonances from the 
regions of motion with a "running" slow angle ip. 

We shall distinguish the relevant saddles as the saddles 
with the lower action/energy (using the subscript "sZ") 
and upper action/energy (using the subscript "sit"). The 
positions of the saddles in the I — ip plane are defined by 
the following equations (cf. Figs. 6 and 7): 

g = sign(q n (I suM )) = (-1)1*1, (16) 
j) sl =n(l + g)/2, $ m = ir(l-g)/2, 
X g (I sl )=X_ g (I su )=0, 



dX g (I s i) 



>o, 



dX_g{I SU ) 



<o, 



dl s i dl su 
where X±(I) are defined in Eq. (15) while I s i and I s . 



are closer to I{E m ) than any other solution of (16) (if 
any) from below and from above, respectively. 

Given that the values of h relevant to the minima of the 
spikes are small in the asymptotic limit $ — > 0, one may 
neglect the last term in the definition of X T in Eq. (15) 
in the lowest-order approximation, so that the equations 
X T = reduce to the simple resonance condition 



nu>(I s 



(17) 



Substituting here Eq. (7) for u>, we obtain the explicit 
expressions for the energies in the saddles: 



(18) 



AE™ ee J$ 2 -64 



exp 



I , uif < nuj m . 



The corresponding actions I su ,si are expressed via E sus i 
by means of Eq. (13). 

For ujf « nw m , the values of E su ,si (18) lie in the range 
where the expression (12) for q„ docs hold true. This will 
be explicitly confirmed by the results of the calculations 
based on this assumption. 

Using (16) for the angles and (18) for the energies, 
and the asymptotic expressions (7), (12) and (13) for 
to(E), q n (E) and 1(E) respectively, and allowing for the 
resonance condition (17), we obtain explicit expressions 
for the values of the Hamiltonian (4) in the saddles: 



H, 



-Hsu — 



7T 



2A£W-$ln( ¥ -^ Ty ) +W2. 



(19) 

As the analysis of simulations suggests (see the item 1 
in the end of Sec. II) and as it is rigorously shown in the 
next subsection, one of the main conditions which should 
be satisfied in the spikes is the overlap in phase space be- 
tween the separatrices of the nonlinear resonances, called 
separatrix reconnection 0,0,0- Given that the Hamil- 
tonian H is constant along any trajectory of the system 
(4), the values of H in the lower and upper saddles of the 
reconnected separatrices are equal to each other: 



H s i — H su , 



(20) 



that may be considered as the necessary and sufficient 
|30| condition for the reconnection. Taking into account 
that H s i = —H su (see (19)), it follows from (20) that 



Hsi = Hsu = 0. (21) 
Explicitly, the relations in (21) reduce to 



8 



h(u f ) 



V2n 



$ In 



A£« = ,/$ 2 _64exp(- — ). 
V u f 

IT 

< cj m — ujf jn <C cj 7) 



(22) 



21n(8/$) : 



ri = 1, 3, 5, 



The function h(wf) (22) monotonously decreases to zero 
as u>f grows from to nu m , where the line abruptly stops. 
Fig. 10 shows the portions of the lines (22) relevant to 
the left wings of the 1st and 2nd spikes (for $ = 0.2). 



C. Barrier chaotic layers 

The next step is to find a minimal value of h for which 
the resonance separatrix overlaps the chaotic layer re- 
lated to a potential barrier. With this aim, we study how 
the relevant outer boundary of the chaotic layer behaves 
as h and u> / vary. Assume that the relevant w/ is close to 
nui rn while the relevant h is sufficiently large for lo(E) to 
be close to w m at all points of the outer boundary of the 
layer (the results will confirm these assumptions). Then 
the motion along the regular trajectory infinitesimally 
close to the layer boundary may be described within the 
resonance approximation (4). Hence the boundary may 
also be described as a trajectory of the resonant Hamilto- 
nian (4). This is explicitly proved in the Appendix, using 
the separatrix map analysis that allows for the validity of 
the relation cu(E) w uj m for all E relevant to the bound- 
ary of the chaotic layer. The main results are presented 
below. For the sake of clarity, we present them for each 
layer separately, although they are similar in practice. 



1. Lower layer 

Let u>f be close to any of the spikes minima. 

One of the key roles in the formation of the upper 
boundary of the layer is played by the angle-dependent 
quantity (^|sin(V>)| which we call the generalized sepa- 
ratrix split (GSS) for the lower layer, alluding to the 
conventional separatrix split 0] for the lower layer Si = 
\e^ low \cjf)\h with e^ ™) given by Eq. (All) 31}. Ac- 
cordingly, we use the term "lower GSS curve" for the 
following curve in the I — ip plane: 



$) = j(i5W + 6l \ s in$)\). 



a. Relatively small h 



(23) 



If h < hir(ujf), where the critical value hcr{wt) is 
determined by Eq. (41) (its origin will be explained fur- 
ther), then there are differences in the boundary forma- 



tion for the frequency ranges of odd and even spikes. We 
describe these ranges separately. 

a.l. Odd spikes 

In this case, the boundary is formed by the trajectory 
of the Hamiltonian (4) tangent to the GSS curve (see 
Fig. 16(a); cf. also Figs. 8(a), 9(b), 9(c)). There are two 
tangencies in the angle range ] — 7r, tt[: they occur at the 
angles ±V>^ where V>j is determined by Eq. (A21). 

In the ranges of h and a->/ relevant to the spike mini- 
mum, the asymptotic expressions for Si and ijyf 1 are: 



Si = V2nh, 



(0 



(-1)1 



81n(l/$) 



.1 -(-!)[' 



(24) 
(25) 



Hence, the asymptotic value for the deviation of the tan- 
gency energy from the lower barrier reduces to: 

IT 3 ' 2 h 



E, 



E, 



(i) 



Si sin(^ 



2 v/ln(l/$)/n 



(26) 



The minimal energy on the boundary, -E^L, corre- 
sponds to ip = or 7r for even or odd values of [n/A] 
respectively. Thus, it can be found from the equality 



H 



(l(^L)^ = 7r(l-(-l)W)/l 



HL 



r(0 



(27) 



At $ — * 0, Eq. (|2"7| yields the following expression for 
the minimal deviation of energy on the boundary from 
the barrier: 



r(0 = 



,(/) 

'min 
^3/2 



E y ' - E 



{1) = {E^-E^)I^- 
h 



2y/e ^/ln(l/$) jn 



(28) 



In the context of global chaos onset, the most impor- 
tant property of the boundary is that the maximal devia- 
tion of its energy from the barrier, Smax, greatly exceeds 



both and Sr. As h 



the maximum of the 



boundary approaches the saddle u sV . 

a. 2. Even spikes 

In this case, the Hamiltonian (4) possesses saddles "s" 
in the close vicinity to the lower barrier (see Fig. 16(b)). 
Their angles differ by it from those of "sf : 



.!-(-!)[* 



27rrn, 



(29) 



m = 0,±l,±2,. 



while the deviation of their energies from the barrier still 
lies in the relevant (resonant) range and reads, in the 
lowest-order approximation, 



2V2ln(MW)' 



(30) 



The lower whiskers of the separatrix generated by these 
saddles intersect the GSS curve while the upper whiskers 
in the asymptotic limit do not intersect it (Fig. 16(b)). 
Thus, it is the upper whiskers of the separatrix which 
form the boundary of the chaotic layer in the asymptotic 
limit. The energy on the boundary takes the minimal 
value right on the saddle "s" , so that 




-3jl -2n 



7T 



h 



5 {l) = 5 = 
mm S 2v^ln(ln(l/$))" 



(31) 



Similar to the case of the odd spikes, the maximal 
(along the boundary) deviation of the energy from the 



maximum of the boundary approaches the saddle "sf 



barrier greatly exceeds both <5^ ; in and Si. As h — > her , the 



b. Relatively large h 

If h > /icr(w/), the previously described trajectory 
(the tangent one or the separatrix, for the odd or even 
spike ranges respectively) is encompassed by the separa- 
trix of the lower nonlinear resonance and typically forms 
the boundary of the major stability island inside the 
lower layer (reproduced periodically in ip with the period 
27r). The upper outer boundary of the layer is formed by 
the upper part of the resonance separatrix. This may be 
interpreted as the absorption of the lower resonance by 
the lower chaotic layer. 



2. Upper layer 

Let u>f be close to any of the spikes minima. 

One of the key roles in the formation of the lower 
boundary of the layer is played by the angle-dependent 
quantity 8 U \ sm(ip)\ which we call the generalized separa- 
trix split (GSS) for the upper layer; S u is the separatrix 
split for the upper layer: 5 U = \e( up \u) f)\h with e^ upS> 
given by Eq. (A43). Accordingly, we use the term "up- 
per GSS curve" for the following curve in the I — tp plane: 



J^ s (V0 = /(i?r-^sin(V0|). 



(32) 



a. Relatively small h 
If ft, < hi"\ujf), where the critical value hi/(u>f) 



is 



determined by Eq. (42) (its origin will be explained fur- 
ther), then there are some differences in the boundary 
formation in the frequency ranges of odd and even spikes: 
for odd spikes, the formation is similar to the one for even 
spikes in the lower-layer case and vice versa. 




FIG. 8: (a). Chaotic layers (shaded in green and blue, for the 
upper and lower layers respectively) in the plane of action / 
and slow angle ip, as described by our theory. Parameters 
are the same as in Figs. 5(b) and 6(b). The lower and up- 
per boundaries of the figure box coincide with I(E^ ) and 
I(E^) respectively. The resonance separatrices are drawn 
by the cyan and red solid lines (for the lower and upper res- 
onances respectively). Dashed green and blue lines mark 
the curves I = 7^ s (^) = I(E = E ( b 1} + ft|sin(^)|) and 



/ = Icssi^) = = ~ <5u| sin(i/>)|) respectively, where 
Si and S u are the values of the separatrix split related to the 
lower and upper barrier respectively. The upper boundary of 
the lower layer is formed by the trajectory of the resonant 
Hamiltonian system (4) tangent to the curve I = Jqs S (^). 
The lower boundary of the upper layer is formed by the lower 
part of the upper (red) resonance separatrix. The periodic 
closed loops (solid blue lines) are the trajectories of the sys- 
tem (4) tangent to the curve I ossify: they form the bound- 
aries of the major stability islands inside the upper chaotic 
layer, (b). Comparison of the chaotic layers obtained from 
computer simulations (dots) with the theoretically calculated 
boundaries (solid lines) shown in the box (a). 



In this case, the Hamiltonian (4) possesses saddles "s" 
in the close vicinity to the upper barrier, analogous to 
the saddles "s" near the lower barrier in the case of even 
spikes. Their angles are shifted by it from those of "s" : 



Ips = 1p s + 7T = 7T- 

m = 0,±l,±2,.. 



1 



(-!)[; 



+ 27TTO, (33) 



a. 1. Odd spikes 



The deviation of their energies from the upper barrier 
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coincides, in the lowest-order approximation, with 6 S : from the barrier: 



Sz = S s 



2V2ln(ln(l/$))' 



(34) 



The upper whiskers of the separatrix generated by 
these saddles intersect the upper GSS curve while the 
lower whiskers in the asymptotic limit do not intersect 
it. Thus, it is the lower whiskers of the separatrix which 
form the boundary of the chaotic layer in the asymptotic 
limit. The deviation of energy from the upper barrier 
takes its minimal (along the boundary) value right on 
the saddle "s" , 



c(«) 



Ss 



2V2ln(ln(l/$))' 



(35) 



The maximal (along the boundary) deviation of the 

energy from the barrier greatly exceeds both 5^ n and S u . 

As h — > , the maximum of the boundary approaches 
the saddle "sm" . 

a. 2. Even spikes 

The boundary is formed by the trajectory of the Hamil- 
tonian (4) tangent to the GSS curve. There are two tan- 
gencies in the angle range ] — ir, n[: they occur at the 



angles ±^"' ) where ij)\ a> is determined by Eq. (A41). 

In the ranges of h and uif relevant to the spike mini- 
mum, the expressions for 8 U and •ip^ in the asymptotic 
limit $ — > are similar to the analogous quantities in the 
lower-layer case: 



-,(«) 



5 U = V2irh, 
$"> = -(-!)[*] 



(36) 



^ +7 t 1 + ( : 1)[TJ .(37) 



s Mi) 



Hence, the asymptotic value for the deviation of the tan- 
gency energy E[ u) from the upper barrier reduces to: 



^3/2 



2 v/ln(l/$)/n 



(38) 



The maximal energy on the boundary, £ma X , corre- 
sponds to i> = tt(1 + (-l)["/ 4 l)/2. Thus, it can be found 
from the equality 

H(I = I(E<&),$ = tt(1 + (-l)["/ 4 l)/2) = 

= H(I^=I(E^),^). (39) 

At <I> — > 0, Eq. (39) yields the following expression for 
the minimal (along the boundary) deviation of energy 



? («) = E (2) 



b 

^3/2 



-E^ = (E^-E^)/V-e = 



2eV2 ^/EJTJWj I r 



(40) 



b. Relatively large h 



If h > hir\tj->f) (cf. Fig. 8(a)), the previously de- 
scribed trajectory (tangent one or the separatrix, for the 
even and odd spikes ranges respectively) is encompassed 
by the separatrix of the upper nonlinear resonance and 
typically forms the boundary of the major stability island 
inside the upper layer (reproduced periodically in ifj with 
the period 2tt). The lower outer boundary of the layer 
is formed in this case by the lower part of the resonance 
separatrix. This may be interpreted as the absorption of 
the upper resonance by the upper chaotic layer. 

The description of chaotic layers given above and, in 
more details, in the Appendix is the first main result 
of this paper. It provides a rigorous base for our intu- 
itive assumption that the minimal value of h at which the 
layers overlap corresponds to the reconnection of the non- 
linear resonances with each other while the reconnected 
resonances touch one of the layers and touch/overlap an- 
other layer. It is remarkable also that we have managed 
to obtain the quantitative theoretical description of the 
chaotic layers boundaries in the phase space, including 
even the major stability islands, that well fits the results 
of simulations (see Fig. 8(b)). 



D. Onset of global chaos: the spikes minima 

The condition for the merger of the lower resonance 
and the lower chaotic layer may be written as 



H(I = I(E = E^+5^ n ),i; = 0)=H sl 



(0 



(41) 



The condition for the merger of the upper resonance 
and the upper chaotic layer may be written as 



H(I = I(E = El 



(2) 



), i> = tt) = H su . 



(42) 



For the global chaos onset related to the spike mini- 
mum, either of Eqs. (41) and (42) should be combined 
with the condition of the separatrix reconnection (20). 
Let us seek first only the leading terms of h s = h s (§) 
and uj s = Ld s (&)- Then (20) may be replaced by its 
lowest-order approximation (21) or, equivalently, (22). 
Using also the lowest-order approximation for the barri- 
ers {E { ^' 2) w 1/2 we reduce Eqs. (41), (42) respec- 
tively to 



H(I 
H(I 



I(E 
I(E 



1/2- $ 
1/2 + $ 



s(0 

min 



),V> 



0) 

TT) 



0, (43) 
= 0, (44) 
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where £® n is given by (28) or (31) for the odd or even 

spikes respectively while <5^ in is given by (35) or (40) for 
the odd or even spikes respectively. 



The solution (h. 



(0 ,,(Q 



of the system of equations 

,(«) , 



(22), (43) and the solution (h s w ', ujf') of the system of 
equations (22), (44) turn out identical to the leading or- 
der. For the sake of brevity, we derive below just 



(h 



(0 ,,W 



), denoting the latter, in short, as (h s ,uj s ) [3c 



The system of algebraic equations (22) and (43) is still 
too complicated to find its exact solution. However, we 
need only the lowest- order solution - and this simplifies 
the problem. Still, even this simplified problem is not 
trivial, both because the function h s (<&) turns out to be 
non-analytic and because AE^ 1 ' in (22) is very sensitive 
to ujf in the relevant range. Despite these difficulties, we 
have found the solution in a self- consistent way, as briefly 
described below. 

Assume that the asymptotic dependence /i s (<£>) is: 



ln(4e/$) 



(45) 



where the constant a may be found from the asymptotic 
solution of (22), (43), (45). 

Substituting the energies E = 1/2 — $ 



1/2 + $ - stL in (7) and taking into account (28), (31), 
(35), (40) and (45), we find that, both for the odd and 
even spikes, the inequality 



€1 and£ 



uj m - u)(E) -C uj m (46) 
holds in the whole relevant range of energies, i.e. for 



Substituting (48) into the expression for AE^ in (22), 
using (45) and keeping only the leading terms, we obtain 



AE W = -4e c " 2 , 



2^2 
c = a. 



(49) 



Substituting AE^ from (49) into Eq. (22) for h(u f ) 
and allowing for (45) once again, we arrive at the tran- 
scendental equation for c: 



In 



1 + X(c) 



2x(c) 



,1-X(c) 

X(c) ee a/1 - 4e c " 2 . 

The approximate numerical solution of Eq. (50) is: 



(50) 



0.179. 



(51) 



Thus, the final leading-order asymptotic formulas for 
ojf and h in the minima of the spikes are the following: 



Uso = w; = n 



hso 



21n(f)' 
c $ 



(52) 



'V21h(t)' 
71 = 1,3,5,..., $^0, 

where the constant c ~ 0.179 is the solution of Eq. (50). 

The rigorous derivation of the explicit asymptotic for- 
mulas for the minima of h gc (ojf) is the second main 
result of this paper. These formulas allow one to im- 
mediately predict the parameters for the weakest pertur- 
bation which may lead to global chaos. 



AE e + 5 {l) , $ - 5 {u) } 



(47) 



Thus, the resonant approximation is valid in the whole 
range (47). Eq. (12) for q n (E) is valid in the whole 
relevant range of energies too. 

Consider Eq. (43) in a more explicit form. Namely, we 
express Vf from (43), using Eqs. (4), (12), and (13), and 
using also (28)/(31) for odd/even spikes, and (45): 



21n (f) 



hV2 



O 



hT(4e/$), 



(48) 



We emphasize that the value of 5)J iri enters explicitly only 
the term 0{. . .) while, as it is clear from the considera- 
tion below, this term does not affect the leading terms in 
(h s , Uj). Thus, S[^ in does not affect the leading term of u> s 
at all, while it affects the leading term of h s only implic- 
itly: <5£. n lies in the range of energies where nq n (E) ps y/2. 
This latter quantity is present in the second term in the 
curly brackets in (48) and, as it is clear from the further 
consideration, h s is proportional to it. 



E. Numerical example and next-order corrections 

For $ = 0.2, the numerical simulations give the follow- 
ing values for the frequencies in the minima of the first 
two spikes (see Fig. 4): 



to 



(i) 



0.4005 ± 0.0005, 



,(2) 



1.24 ±0.005. (53) 



The values by the lowest-order formula (52) are: 



X 1 ) 



lj s0 « 0.393, 



,( 2 ) 



1.18, 



(54) 



in rather good agreement with the simulations. 

The next-order correction for uj s can be immediately 
found from Eq. (48) for uif and Eq. (52) for h S Q, so that 



717T 1 



n = 1,3, 5, ... 



0.09 



21n(t) 



21n(f) 



,(55) 
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The formula (55) agrees with the simulations even bet- 
ter than the lowest-order approximation: 



J s \ } ps 0.402, ps 1.21. 



(56) 



For h in the spikes minima, the simulations give the 
following [34[ values (see Fig. 4): 



h^ ps 0.0049, re 0.03. (57) 

The values by the lowest-order formula (52) are: 



0.0032, 



The theoretical value h^J gives a reasonable estimate 



0.01. 



(58) 



for the simulation value hi 1 ^. The theoretical value h 



,(2) 
sO 



gives the correct order of magnitude for the simulation 

( 2) 

value hg ■ Thus, the accuracy of the lowest-order for- 
mula (52) for h a is much lower than that for u> s : this is 
due to the steepness of h gc (u)f) in the ranges of spikes 
(the steepness, in turn, is due to the flatness of the func- 
tion ui(E) near its maximum). Moreover, as the number 

of the spike j increases, the accuracy of the lowest-order 

(i) 

value h sQ significantly decreases. The latter can be ex- 
plained as follows. For the next-order correction to h^f, 
the dependence on $ reads as: 



a sl 



0) 



,0) 

l s0 



oc 



h (J) ~ln(4e/<I>)' 



(59) 



At least some of the terms of this correction are posi- 
tive and proportional to h^Q (e.g. due to the difference 
between the exact equation (15) and its approximate ver- 
sion (17)) while is proportional to n = 2j — 1. Thus, 
for $ = 0.2, the relative correction for the 1st spike is 
~ 0.25 while the correction for the 2nd spike is a few 
times larger i.e. ~ 1. But the latter means that, for 
$ = 0.2, the asymptotic theory for the 2nd spike cannot 

(2) 

pretend to be a quantitative description of h s , but only 
provide the correct order of magnitude. Besides, if n > 1 
while $ exceeds some critical value, then the search of 
the minimum involves Eq. (66) rather than Eq. (20), 
as explained below in Sec. IV (cf. Figs. 10(b) and 11). 
Altogether, this explains why hs is larger than h s $ only 

by 50% while is larger than h% } by 200%. 

The consistent explicit derivation of the correction to 
h^Q is complicated. A reasonable alternative may be 
a proper numerical solution of the algebraic system of 
Eqs. (20) [H and (41) for the odd spikes or (42) for 

the even spikes [H & To this end > in E q s - ( 20 ) HI 
and (41)/(42) we use: (i) the exact values of the saddle 
energies obtained from the exact relations (16) instead of 
the approximate relations (17); (ii) the exact formulas (5) 



and (6) for ui{E) instead of the asymptotic expression (7); 
(iii) the exact functions q n (E) instead of the asymptotic 
formula (12); (iv) the relation (27) and the calculation 

of the "tangent" state ($ ,4°) by Eqs. (All), (A22) 
for the odd spikes, or relation (39) and the calculation 

of the "tangent" state (tjj { t u) , I { t u) ) by Eqs. (A41)-(A43) 
for the even spikes. Note that, to find the exact function 
q n (E), we substitute into the definition of q n {E) in (4) 
the explicit [38[ solution for q(E,ip): 



q(E, ip) — arcsin 



V 



1—7/ 

q(E,ip) = 7r- q(E,n- i/j) 
q(E,i>)=q(E,27r-iP) 

T] = i(V2^-$ + l)sn 2 (—ill 

2 V 7T 





7T 


for ij) e 


M 



7T 

2' 77 



for ip £ 
for ip £ [it, 2tt] , 

(60) 



where sn(x) is the elliptic sine [25j with the same modulus 
k as the full elliptic integral K defined in (5), (6). 
The numerical solution described above gives: 



(oj^A sa 0.401, (h^) « 0.005, 

V / num \ / num. 

("'■') 

\ / nu 



(61) 



1.24 



h? 



0.052 



The agreement with the simulation results is: (i) ex- 
cellent for uo s for the both spikes and for h s for the 1st 
spike, (ii) reasonable for h s for the 2nd spike. Thus, if 
is moderately small, a much more accurate prediction for 
h s than that by the lowest-order formula is provided by 
the numerical procedure described above. 



IV. THEORY OF THE SPIKES WINGS 

The goal of this section is to find mechanisms responsi- 
ble for the formation of the spikes wings (i.e. the function 
h gc (u)f) in the ranges of u>f slightly deviating from u>i ) 
and to provide for their theoretical description. 

Before developing the theory, we briefly analyze the 
simulation data (Fig. 4), concentrating on the 1st spike. 
The left wing of the spike is smooth and nearly straight. 
The initial part of the right wing is also nearly straight 
[39l | though less steep. But, at some small distance from 
cui 1 ^ , its slope changes] ump- wise by a few times: com- 
pare the derivative [39|] dh gc /duif ps 0.1 at ut slightly 

exceeding wj 1 ' ~ 0.4 (see the left inset in Fig. 4) and 
dhg C /dujf ps 0.4 at u)f = 0.45 -j- 0.55 (see the main part 
of Fig. 4). Thus, even prior to the theoretical analysis, 
one may assume that there are a few different important 
mechanisms responsible for the formation of the wings. 

Consider the arbitrary jth spike. We have shown in 
the previous section that, in the asymptotic limit $ — > 0, 
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FIG. 9: Illustrations of the mechanisms of the formation of 
the 1st spike wings and of the corresponding theoretical lines 
in Fig. 10(a). Boxes (a), (b) and (c) illustrate the lines 
by Eqs. (20), (41) and (64) respectively: the correspond- 
ing perturbation parameters are (luj = 0.39, h = 0.0077), 
{cjf = 0.41, h = 0.00598) and (u f = 0.43, h = 0.01009) re- 
spectively. Resonance separatrices are drawn by red and cyan. 
The dashed lines show the functions Iq S q('4>) and ^gss(^). 
The black line in (c) is the trajectory of the resonant Hamil- 
tonian system (4) which is tangent to both dashed lines. 



the minimum of the spike corresponds to the intersec- 
tion between the lines (20) and (41)/(42) for odd/even 
spikes. We recall that: (i) Eq. (20) corresponds to the 
overlap in phase space between nonlinear resonances of 
the same order n = 2j — 1; (ii) Eq. (41)/(42) corre- 
sponds to the onset of the overlap between the resonance 
separatrix associated with the lower/upper saddle and 
the chaotic layer associated with the lower/upper poten- 
tial barrier; (iii) for = u>g , the condition (41)/(42) 
guarantees also the overlap between the upper/lower res- 
onance separatrix and the chaotic layer associated with 
the upper/lower barrier [33| . 

If Uf becomes slightly smaller than loI the resonances 
shift closer to the barriers while moving apart from each 



other. Hence, as h increases, the overlap of the reso- 
nances with the chaotic layers associated with the bar- 
riers occurs earlier than with each other. Therefore, at 
< Ws — <C u) m , the function h gc (ujf) should approx- 
imately correspond to the reconnection of resonances of 
the order n = 2j — l (Fig. 9(a)). Fig. 10(a) demonstrates 
that even the asymptotic formula (22) for the separatrix 
reconnection line fits the left wing of the 1st spike quite 
well while the numerically calculated line (20) agrees with 
the simulations perfectly. 

If uif becomes slightly larger than uj s then, on the 
contrary, the resonances shift closer to each other and 
more far from the barriers. Therefore, the overlap of 
resonances with each other occurs at smaller h than the 
overlap between any of them and the chaotic layer asso- 
ciated with the lower/upper barrier (cf. Figs. 5(c) and 
5(d) as well as 6(c) and 6(d)). Hence, it is the latter over- 
lap which determines the function h gc (ujf) in the relevant 
range of w/ (Fig. 9(b)). Fig. 10 shows that h gc (ujf) is in- 
deed well approximated in the close vicinity to the right 
from Wj by the numerical solution of Eq. (41)/(42) for 
an odd/even spike and, for the 1st spike and the given 
<&, even by its asymptotic form, 



h = h{w f ) = (62) 
${21n(f)+l n (|±A|^)}-2A^ 



2V2 



A£ (1) = ,/$2-64exp( ), 



n = 2j-l, 



The mechanism described above determines h„ 
only in the close vicinity of Us . 



If ujf/n becomes too 
close to io m or exceeds it, then the resonances are not of 
immediate relevance: they may even disappear or, if they 
still exist, their closed loops shrink, so that they cannot 
anymore provide for the connection of the chaotic lay- 
ers in the relevant range of h. At the same time, the 
closeness of the frequency to uj m still may give rise to 
a large variation of action along the trajectory of the 
Hamiltonian system (4). For the odd/even spikes, the 
boundaries of the chaotic layers in the asymptotic limit 
$ — * are formed in this case by the trajectory of (4) 
which is tangent to the lower/upper GSS curves (for the 
lower/upper layer) or by the lower/upper part of the sep- 
aratrix of (4) generated by the saddle "s"/"s" (for the 
upper/lower layer). Obviously, the overlap of the layers 
occurs when these trajectories coincide with each other, 
that may be formulated as the equality of H in the cor- 
responding tangency and saddle: 



r(«) „7,(«) 



for j = 1,3,5,. 



H(I s ^ s )=H(ir,i>r>) ^ j = 2, 4, 6, 



(63) 
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FIG. 10: The 1st (a) and 2nd (b) spike in h gc (u)f): com- 
parison between the results of the numerical simulations (the 
lower boundary of the shaded area) and the theoretical es- 
timates. The estimates are indicated by the corresponding 
equation numbers and are drawn by different types of lines, 
in particular the dashed lines represent the explicit asymptote 
for the solid line of the same color. 



Note however that, for moderately small $, the tangencies 
may be relevant both to the lower layer and to the upper 
one (see the Appendix). Indeed, such a case occurs for 
our example with <t> = 0.2: see Fig. 9(c). Therefore, the 
overlap of the layers corresponds to the equality of H in 
the tangencies: 



H(ll l \^) = H(lt\$ u) ). 



(64) 



To the lowest order, Eq. (63) and Eq. (64) read as: 



h = h{u f ) = 



\/2$ In I 



4c \ 
$ ) 



21n(f) 



(65) 



Both the line (64) and the asymptotic line (65) well 
agree with the part of the right wing of the 1st spike 
situated beyond the immediate vicinity of the minimum 
from the right side, namely, to the right from the fold 
at ojf ~ 0.42 (Fig. 10(a)). The fold corresponds to the 
change of the mechanisms of the chaotic layers overlap. 
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FIG. 11: Amplitude of the 3rd Fourier harmonic as a function 
of action (solid red line). The dashed black line shows the zero 
level. Its intersection with the solid red line is marked by the 
circle. The green line indicates the value of action where 
§3 = 0. The inset illustrates the line (66) in Fig. 10(b): the 
GSS curve touches the horizontal line / — /q 3 — o- 



If $ is moderately small while n > 1, the description 
of the far wings by the numerical lines (20) and (64) may 
be still quite good but the asymptotic lines (22) and (65) 
cannot pretend to describe the wings quantitatively any- 
more (Fig. 10(b)). As for the very minimum of the spike 
and the wings in the close vicinity to it, one more mecha- 
nism may become relevant for their formation in this case 
(Figs. 10(b) and 11). This mechanism may be explained 
as follows. If n > 1, then q n (E) becomes zero in the close 
vicinity (~ $ ) of the relevant barrier (the upper/lower 
barrier, in the case of even/odd spikes: cf. Fig. 11). As 
follows from the equations of motion (14), no trajectory 
can cross the line I — Iq n =o- In the asymptotic limit 
<f> — > 0, provided h is from the relevant range, almost the 
whole GSS curve is farer from the barrier than the line 



/ = Iq n =o, and the latter becomes irrelevant. But, for a 
moderately small the line may separate the whole GSS 
curve from the rest of the phase space. Then the reso- 
nance separatrix cannot connect to the GSS curve even 
if there is a state on the latter curve with the same value 
of H as on the resonance separatrix. For a given u)f, the 
connection requires then a higher value of h: for such a 
value, the GSS curve itself crosses the line / = I qn= o- In 
the relevant range of h, the resonance separatrix passes 
very close to this line, so that the connection is well ap- 
proximated by the condition that the GSS curve touches 
this line (see the inset in Fig. 11): 



5 U = - E q2] _ 1=0 for j = 2, 4, 6, 



(66) 



This mechanism is relevant for the formation of the mini- 
mum of the 2nd spike at $ = 0.2, and in the close vicinity 
of the spike, on the left (Fig. 10(b)). 

Finally, let us explicitly find the universal asymptotic 
shape of the spike in the vicinity of its minimum. 

First, we note that the lowest-order expression (62) 
for the spike between the minimum and the fold can be 
written as the half-sum of the expressions (22) and (65) 
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(which represent the lowest-order approximations for the 
spike to the left of the minimum, and to the right of the 
fold respectively). Thus, all three lines (22), (62) and 
(65) intersect in one point. The latter means that, in the 
asymptotic limit $ — ► 0, the fold merges with the mini- 
mum: Uf and h in the fold asymptotically approach u> s 
and h s respectively. Thus, though the fold is a generic 
feature of the spikes, it is not of main significance: the 
spike is formed basically by two straight lines. The ra- 
tio between their slopes is universal. So, introducing a 
proper scaling, we reduce the spike shape to a universal 
function (Fig. 12): 

h(Au)f) = h (lw) (Au f ) = 1 - y/l - Ae c ~ 2 ACo s ps 

w 1 - 0.593Aw/ for Au> f < 0, 
h(Aui f ) = U rw \ACb s ) = 1 + Aw f for Aui f > 0, 

(67) 

~ h (fol d){AQf) = hW{te f ) + hl™\*u, f ) _ 

1 - a/1 — 4c c ~ 2 
= 1 + —- Au f ps 1 + 0.203Aw/, 

h = —, Au f = ^ , $^0, 

where u s o and h s o are the lowest-order expressions (52) 
respectively for the frequency and amplitude in the spike 
minimum, u> s i is the expression (55) for the frequency in 
the spike minimum, including the first-order correction, 
and c is the constant (51). 

Beside the left and right wings of the universal shape 
(the solid lines in Fig. 12), we also present in (67) the 
function U iold \Auj f ) (the dashed line in Fig. 12): it 
purposes to show that, on one hand, the fold asymptoti- 
cally merges with the minimum but, on the other hand, 
the fold is generic and the slope of the spike between the 
minimum and the fold has a universal ratio to any of the 
slopes of the major wings. 

Even for a moderately small <&, like in our example, 
the ratios between the three slopes related to the 1st 
spike in the simulations are reasonably well reproduced 
by those in Eq. (67): cf. Figs. 10(a) and 12. It follows 
from (67) that the asymptotic scaled shape is universal 
i.e. independent of <I>, n or any other parameter. 

The description of the wings of the spikes near their 
minima, in particular the derivation of the spike universal 
shape, is the third main result of this paper. 

V. GENERALIZATIONS AND APPLICATIONS 

The new approach for the treatment of separatrix chaos 
opens a broad variety of important generalizations and 
applications, some of which are discussed below. 

1. It may be applied to any separatrix layer, includ- 
ing in particular the conventional single- separatrix 




FIG. 12: The universal shape of the spike minimum (67) 
(solid lines). The dashed line indicates the universal slope 
of the spike in between the minimum and the fold, which 
have merged in the universal (asymptotic) function (67). 

case. This is possible due to the characteristic de- 
pendence of the frequency of eigenoscillation on en- 
ergy in the vicinity of any separatrix: the frequency 
keeps nearly a constant value even if the deviation 
of the energy from the separatrix strongly varies 
within a given scale of the deviation. To the best 
of our knowledge, the latter idea was not exploited 
before, which is why probably it was thought to 
be impossible to analyze the dynamics of the sepa- 
ratrix map in an explicit form. There were var- 
ious estimates of the layer width in energy (see 
EL S IS EL IE H2] and references therein) but the 
quantitative analysis of the layer boundaries in the 
phase space was never done in the non-adiabatic 
case. In contrast, our approach does analyze the 
dynamics of the separatrix map, and this allows 
us in particular to quantitatively describe the layer 
boundaries in case when the perturbation is reso- 
nant with the eigenoscillation in the relevant range 
of energy. It follows from such a description that, 
for a given small amplitude h of the perturba- 
tion, the maximal layer width in energy is much 
larger than it is assumed by former theories (cf. 
[3)13) B EL where the maximal width is as- 
sumed to be ~ h). Thus, we can quantitatively 
describe large jumps and peaks of the layer width 
as a function of the perturbation amplitude and fre- 
quency respectively [EES]- Moreover, a rough es- 
timate on the basis of our approach indicates that, 
for many classes of systems, the relative range of 
such jumps/peaks (i.e. the ratio between the up- 
per and lower levels of the jump/peak) diverges in 
the asymptotic limit h — > 0. 

2. Apart from the description of the boundaries, the 
approach allows us to describe the chaotic transport 
within the layer. In particular, it may allow us 
to calculate a positive Lyapunov exponent and to 
describe diffusion. 

3. Our approach may be generalized for a non- 
resonant perturbation. The resonance approxima- 
tion is not valid then but there still remains the 
property of a nearly constancy of the frequency of 
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eigenoscillation within an arbitrary given scale of 
the deviation of energy from the separatrix. This 
property may allow us to explicitly describe the dy- 
namics of the separatrix map for any frequency of 
perturbation which is less than or of the order of 
the resonance frequency. 

4. Apart from Hamiltonian systems of the one and 
a half degrees of freedom and corresponding Za- 
slavsky separatrix maps, our approach may be use- 
ful in the treatment of other chaotic systems and 
separatrix maps (see [32| for the most recent ma- 
jor review on various types of separatrix maps and 
related continuous chaotic systems). 

As concerns the facilitation of the global chaos onset 
between adjacent separatrices, we mention the following 
generalizations : 

1. The spikes in h gc {ujf) may occur for an arbitrary 
Hamiltonian system with two or more separatrices. The 
asymptotic theory can be generalized accordingly. 

2. The absence of pronounced spikes at even harmon- 
ics 2juj m is explained by the symmetry of the poten- 
tial (1): the even Fourier harmonics of the coordinate, 
<72j, are equal to zero. For time-periodic perturbation of 
dipole type, like in Eq. (2), there are no resonances of 
even order due to this symmetry P, H, H, 0, H, Q . If either 
the potential is non- symmetric or the additive perturba- 
tion of the Hamiltonian is not an odd function of the co- 
ordinate, then even-order resonances do exist, resulting 
in the presence of the spikes in h gc (ujf) at uif » 2juj m . 

3. There may be an additional facilitation of global 
chaos onset which is reasonable to call a "secondary" fa- 
cilitation. Let the frequency ujf be close to the frequency 
lu s of the spike minimum while the amplitude h be ~ h s 
but still lower than h gc (uif). Then there are two reso- 
nance separatrices in the I — ip plane which are not con- 
nected by the chaotic transport (cf. Fig. 6(b) and 5(b)). 
This system possesses the zero-dispersion property. The 
trajectories of the resonant Hamiltonian (4) which start 
in between the separatrices oscillate in / (as well as in 
dijj/dt). The frequency Q of such oscillations along a 
given trajectory depends on the corresponding value of 
H analogously to as u> depends on E for the original 
Hamiltonian Hq: u>(H) is equal to zero for the values 
of H corresponding to the separatrices (being equal in 
turn to H s i and H su : see Eq. (19)) while possessing a 
nearly rectangular shape in between, provided the quan- 
tity \H s i — H su \ is much smaller than the variation of H 
within each of the resonances, 

\H s i — H su \ -C H var ~ \H s i — H e i\ ~ \H S u — H eu \, (68) 

where H e i and H eu are the values of H in the elliptic 
point of the lower and upper resonance respectively. The 
maximum of Cj{H) in between H s i and H su is asymptot- 
ically described by the following formula: 



w m ~ — jz ; ; — r- (69) 

In [H var /\H s i — H su \\ 

If we additionally perturb the system in such a way 
that an additional time-periodic term of the frequency 
ujf sa uj m arises in the resonance Hamiltonian, then the 
chaotic layers associated with the resonance separatrices 
may be connected by chaotic transport even for a rather 
small amplitude of the additional perturbation, due to a 
scenario similar to the one described in this paper. 

There may be various types of such additional pertur- 
bation 41]. For example, one may add to H (2) one more 
dipole time-periodic perturbation of mixed frequency (i.e. 
« uj m + Cj m ). Alternatively, one may directly perturb the 
angle of the original perturbation by a low- frequency per- 
turbation, i.e. the time-periodic term in H (2) is replaced 
by the term 

— hqcos(u)ft + A cos(u ft)), (70) 

Ulf Ri L) m , LJf U! m . 

Recent physical problems where a similar situation is 
relevant are: chaotic mixing and transport in a meander- 
ing jet nowj42j and reflection of light rays in a corrugated 
waveguide [43|. 

4. If the time-periodic perturbation is multiplicative 
rather than additive, the resonances become parametric 
(cf. [2l|). Parametric resonance is more complicated and 
much less studied than nonlinear resonance. Neverthe- 
less, the main mechanism for the onset of global chaos 
remains the same, namely the combination of the recon- 
nection between resonances of the same order and of their 
overlap in energy with the chaotic layers associated with 
the barriers. At the same time, the frequencies of the 
main spikes in h gc (u>f) may change (though still being 
related to u) m ). We consider below the example when 
the periodically driven parameter is the parameter $ in 
(1) [44j. The Hamiltonian reads as 

H = p 2 /2+($-sin(g)) 2 /2, (71) 
$ = $ + hcos(idft), $o = const < 1. 

The term ($ — sin(q)) 2 /2 in H (71) may be rewritten as 
($ -sm(q)) 2 /2 + ($o-sin-(q))hcos(ujft)+h 2 cos 2 (ui f t)/2. 
The last term in the latter expression does not affect the 
equations of motion. Thus, we end up with an additive 
perturbation ($o — sm(q))h cos(u> ft). In the asymptotic 
limit <&o — > 0, the nth-order Fourier component of the 
function ($o — sin(q)) can be shown to differ from zero 
only for the orders n — 2,6,10,... Therefore one may 
expect the main spikes in h gc (u>f) to be at frequencies 
twice larger than those for the dipole perturbation (2): 

w£> w 2wC0 a 2 (2j - l)u m , j = 1,2,3, ... (72) 
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FIG. 13: The diagram analogous to that one in Fig. 4 but for 
the system (71) (with $ = 0.2). 



This well agrees with the results of simulations (Fig. 13). 

Moreover, the asymptotic theory for the dipole per- 
turbation may be immediately generalized to the present 
case: it is necessary only to replace the Fourier compo- 
nent of the coordinate q by the Fourier component of the 
function (<£>o — sin(g)): 



($o - sm(q)) n = | ™ 
j = 1,2,3,..., 

(cf. Eq. (12) for q n ). We obtain: 



at n=2(2j-l), 
at ra#2(2j-l), 



(73) 



UJ 



spa 



J sp0 



7T 



"spO — "-spO 

n = 2,6,10,. 



C7T $ 



(74) 



where c is given in Eqs. (50) and (51). 

For $o = 0.2, Eq. (74) gives, for the 1st spike, values 
differing from the simulation data by about 3% for the 
frequency and by about 10% for the amplitude. Thus, the 
lowest-order formulas accurately describe the 1st spike 
even for a moderately small 

5. One more generalization relates to multi- 
dimensional Hamiltonian systems with two or more sad- 
dles with different energies: the perturbation may not 
necessarily be time-periodic, in this case. The detailed 
analysis will be done elsewhere. 

Finally, we point out some analogy between the facili- 
tation of global chaos onset described in our work and the 



so called stochastic percolation in 2D Hamiltonian sys- 
tems described in |45| : the merging of internal and exter- 
nal chaotic zones is also relevant there, like in our case. 
However, both the models and the underlying mecha- 
nisms are very different. Namely, the problem studied 
in [45| | is not of the zero-dispersion type; and the two- 
dimensionality is inherently important. 

Let us turn now to a few detailed examples of applica- 
tions of the facilitation of global chaos onset. 



A. Electron gas in a magnetic superlattice, 
spinning pendulum, cold atoms in an optical lattice 

The first application relates to a classical electron gas 
in a magnetic superlattice [H, Ell E3, El, Ell H3| , where 
the electrons may be considered as non-interacting quasi- 
particles moving on a plane perpendicular to the mag- 
netic field spatially periodic in one of the in-plane di- 
rections (we denote it as x-direction) . Then the electron 
motion in the x-direction is described by the Hamiltonian 
(1) in which q and p are the scaled electron coordinate 
x and generalized momentum p x respectively, while the 
parameter $ is proportional to the reciprocal amplitude 
of the external magnetic field, B -1 , and to the gener- 
alized momentum p„ in the second ( perp endicular to x) 



llperpcni 
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16 L Note that 



in-plane direction: see for details 
p y remains constant during the motion 

If an AC electric field is applied in the x-direction, then 
the model (2)-(l) becomes relevant. The DC- conductivity 
in the cc-direction is proportional to the fraction of elec- 
trons that can take part to the unbounded motion in the 
x-direction. This fraction, in turn, significantly grows as 
the range of energies involved in the unbounded chaotic 
transport increases [l5l.[l6l]. 

If electrons move in vacuum [46| , then it may be possi- 
ble to inject a beam of electrons which possess the same 
velocity. In this case, the parameter <I> in the model (1) 
has some certain value, so that the results obtained in the 
previous sections are directly applicable. The spikes in 
h gc (u)f) mean a drastic increase of the DC-conductivity 
occurring at a very weak amplitude of the AC electric 
field. The frequency of the AC field should be close to 
one of the spike frequencies. The effect is especially pro- 
nounced for the 1st spike i.e. when u>s m 7r/[2 ln(4e/$)]. 

If the electron motion takes place in a semiconductor 
HM Ell E3, EH, then the velocity in the y-direction is 
necessarily statistically distributed. The same concerns 
the parameter <E> then. This might seem to smear the 
effect: cf. [15| where the conventional scenario of the 
onset of global chaos was exploited. However, in the 
case of the "zero-dispersion" scenario suggested in our 
paper, it typically should not be so. Indeed, a statistical 
distribution of the velocity typically decreases exponen- 
tially sharply as the velocity exceeds some characteristic 
value v c : for high temperature T, the Boltzmann dis- 
tribution of energy is relevant and, therefore, v c cx \/T] 
for low temperatures, the Fermi distribution is relevant 
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and, therefore, v c oc \JEp where Ep is the Fermi energy. 
On the other hand, in the range of small <&, the function 
ln(l/$) does not significantly change even if $ changes 
by a few times. Hence, if $ c <C 1 (where $ c denotes the 

$ value corresponding to v c ), then the frequency w^ 1 ' of 
the partial (i.e. for a given value of v y ) lst-order spike 
is nearly the same for most of the velocities in the rel- 
evant range v y ~ v c , and it is approximately equal to 



uj c = 7r/[21n(4e/$ c )]. Similarly, hi^ ~ h c w w c $ c /25. 

Thus, as a function of w/ for fixed h ~ h c , the DC- 
conductivity should have a sharp maximum for wu c . 

If the parameter $ is time-dependent (e.g. if the ex- 
ternal magnetic field has an AC component or/and there 
is an AC electric force perpendicular to the x-direction), 
then the applications may be similar, with the only dif- 
ference that the values of w/ and h in the minima of the 
main spikes differ from those for the additive perturba- 
tion, see Fig. 13 and the related discussion. 

The results of the present paper may also be of direct 
relevance for a pendulum spinning about its vertical axis 
[3], provided the friction is small. The periodic driving 
may be easily introduced mechanically, or electrically if 
the pendulum is electrically charged, or magnetically if 
the pendulum includes a ferromagnetic material. 

Finally, we mention in this subsection that potentials 
similar to U(q) (1), i.e. periodic potentials with two 
different-height barriers per period, may readily be gen- 
erated for cold atoms by means of optical lattices [471 ] . 
The dissipation may be suppressed by means of the de- 
tuning from the atomic resonance 47]. Then the results 
of our paper are also of direct relevance to such systems. 



B. Noise-induced escape 

Consider the noise-induced escape over a potential bar- 
rier in the presence of a non-adiabatic periodic driving. 
For a moderately weak damping, such a driving decreases 
the activation barrier due to the resonant pushing the 
system in the range of resonant energies [48[. If the 
damping is even smaller, the decrease of the activation 
barrier becomes larger, due, however, to a different mech- 
anism, typically related to the chaotic layer associated 
with the separatrix of the unperturbed system 0, l40l ]. 
The lower energy boundary of the layer is smaller than 
the potential barrier energy, so it is sufficient that the 
noise pulls the system right to this boundary (rather than 
to the very top of the potential barrier), after which the 
system may escape over the barrier purely dynamically. 
If the eigenfrequency as a function of the energy pos- 
sesses a local maximum, then the effect may be even more 
pronounced [a, |49(: the decrease of the activation bar- 
rier may become comparable to the potential barrier at 
unusually small amplitudes of the driving, provided the 
driving frequency is close to the extremal eigenfrequency. 
One of the main mechanisms of the latter effect is closely 
related to phenomena discussed in the present paper. In 



the case of escape over two barriers of different heights, 
the effect should become even more pronounced due to 
the mechanism responsible for the spikes of h gc (ujf) stud- 
ied in the present paper. If the potential is periodic, e.g. 
like in optical lattices [47j , the effect may lead to a drastic 
acceleration of the spatial diffusion. 



C. Stochastic web 

Our results may be applied to the stochastic web for- 
mation [H, [5(| HH H2l- If a harmonic oscillator is per- 
turbed by a plane wave whose frequency is equal to the 
oscillator eigenfrequency or its multiple, then the per- 
turbation plays two roles 0, (EH. On one hand, due to 
the resonance with the oscillator, it transforms the struc- 
ture of the phase space of the oscillator, leading to an 
infinite number of cells divided by a unique separatrix. 
It has the form of a web of an infinite radius. On the 
other hand, the perturbation "dresses" this separatrix 
by an exponentially narrow chaotic layer (it is sometimes 
called "stochastic" layer) . Such a web-like layer is called 
stochastic web. It may lead to chaotic transport of the 
system for rather long distances both in coordinate and 
in energy. 

In case when either the resonance is not exact or/and 
the unperturbed oscillator possesses some nonlinearity, 
the perturbation gen erates many separatrices embedded 
into each other (3I l52j] rather than one single infinite web- 
like separatrix. Then a significant chaotic transport in 
energy may arise only if the magnitude of the pertur- 
bation exceeds some critical value corresponding to the 
overlap of chaotic layers associated, at least, with two 
neighbouring separatrices. And, still, the transport in 
energy remains limited since the width of the chaotic 
layer around each separatrix sharpl y d ecreases as the en- 
ergy increases. It can be shown [41j that some types 



of additional time-periodic perturbation lead to a low- 
frequency dipole perturbation of the resonance system 
(cf. the paragraph preceding Eq. (70)). The structure 
of separatrices in the reduced system possesses properties 
similar to that of the system considered in the present pa- 
per. Indeed, in the region between the separatrices, the 
resonance system performs regular oscillations, and the 
frequency of such oscillations, as function of the value 
of the resonance Hamiltonian, is equal to zero at each 
of the separatrices. Thus, it necessarily possesses a lo- 
cal maximum between energies corresponding to any two 
neighbouring separatrices, like in the case considered in 
the present paper. If the additional perturbation has an 
optimal frequency related to one of these local maxima, 
then the overlap of chaotic layers associated with neigh- 
bouring separatrices is greatly facilitated, similarly to the 
case considered in the present paper. Moreover, the lo- 
cal maximum of the eigenfrequency changes from pair to 
pair of separatrices weakly, so that if the magnitude of 
the auxiliary perturbation exceeds the critical value even 
slightly the simultaneous overlap between many chaotic 
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layers may occur. Then, the distance of the chaotic trans- 
port in energy greatly increases. 

Similar applications are relevant for the so called ho- 
mogeneous (sometimes called periodic) stochastic webs 
[H, |5(| and many other web- like stochastic structures [|| . 

Beside classical systems, stochastic webs may arise in 
quantum systems too. It was recently demonstrated, 
both theoretically [53| and experimentally that the 
stochastic web may play a crucial role in quantum elec- 
tron transport in semiconductor superlattices subjected 
to stationary electric and magnetic fields. Due to the spa- 
tial periodicity with a period of about a few nanometers, 
the system possesses narrow minibands in the electron 
spectrum. It turns out that the description of the elec- 
tron transport in the lowest miniband may be approx- 
imated by the model of a classical harmonic oscillator 
driven by a plane wave. The role of the harmonic oscil- 
lator is played by the cyclotron motion while the wave 
arises due to the interplay between the cyclotron motion 
and Bloch oscillations. If the cyclotron and Bloch fre- 
quencies are commensurate, then the phase space of such 
a system is threaded by a stochastic web. This gives rise 
to the derealization of electron orbits, that leads in turn 
to a strong increase of the conductivity However, 
this effect occurs only when the ratio between the elec- 
tric and magnetic fields lies in the exponentially narrow 
regions corresponding to nearly integer ratios between 
the Bloch and cyclotron frequencies. The results of the 
present work suggest a method for a significant increase 
of the width of the relevant regions. If the cyclotron and 
Bloch frequencies are not exactly commensurate, then 
the stochastic web does not arise: rather a set of embed- 
ded separatrices arises provided the effective wave am- 
plitude is sufficiently large. As discussed in the previous 
paragraph, even a rather weak time-periodic driving [55| 
of the optimal frequency may significantly increase the 
area of the phase space involved in the chaotic transport. 
This may provide for an effective control of the electron 
transport in such a system and may be used for develop- 
ing electronic devices that exploit the intrinsic sensitivity 
of chaos (cf. [54(). A similar effect may be used also to 
control transmission through other periodic structures, 
e.g. ultra-cold atoms in optical lattices [56|, [58| and 
photonic crystals (59| . 



VI. CONCLUSIONS 

We have developed a new general approach for the 
treatment of separatrix chaos. This has allowed us to 
create a self-contained theory for the drastic facilitation 
of the global chaos onset between adjacent separatrices 
of a ID Hamiltonian system subject to a time-periodic 
perturbation. Both the new approach and the theory of 
the facilitation are closely interwoven in our paper but, 
at the same time, each of these two items is relevant even 
on its own. That is why we summarize them separately. 

I. The new approach for separatrix chaos. 



The approach is based on the separatrix map analysis 
which uses the characteristic property of the dependence 
of the frequency of eigenoscillation on energy in the vicin- 
ity of a separatrix: the frequency keeps nearly a constant 
value even if the deviation of the energy from the separa- 
trix strongly varies within a given scale of the deviation. 
Due to this, the separatrix map evolves along the major 
part of the chaotic trajectory in regular-like way. The 
deviation of the chaotic trajectory from the separatrix 
may vary along the regular-like parts of the trajectory in 
a much wider range than along the irregular parts. 

In the case of resonant perturbation, we match the 
separatrix map analysis and the resonant Hamiltonian 
approximation. This allows us in particular to find the 
boundaries of the chaotic layers in the phase space, which 
well agrees with computer simulations (Fig. 8). The lat- 
ter theory has been successfully applied by us to the prob- 
lem of the global chaos onset in the double-separatrix 
case, which is summarized in the item II below. Other 
applications and generalizations of the approach include 
in particular the following. 

• It may be applied to the conventional single- 
separatrix case. In particular, our theory predicts 
that the maximal width of the separatrix chaotic 
layer in energy is typically 3> h, in contrast with 
former theories [1, 2, 3, 4, 5] which assume that the 
maximal width is ~ h. 

• It allows to analyze the transport within the layer. 

• It may be generalized for a non-resonant perturba- 
tion and for a higher dimension. 

II. The facilitation of the global chaos onset. 

We have considered in details the characteristic exam- 
ple of a Hamiltonian system possessing two or more sep- 
aratrices, subject to a time-periodic perturbation. The 
frequency lo of oscillation of the unperturbed motion nec- 
essarily possesses a local maximum u) m as a function of 
energy E in the range between the separatrices. It is 
smaller than the frequency ujq of eigenoscillation in the 
stable state of the Hamiltonian system by a factor 

«--(£). -f. p.) 

where At/ is the difference of the separatrices energies, 
while U is the difference between the upper separatrix 
energy and the stable state energy. 

If $ <C 1, the function w(E) is close to uo m for most of 
the energy range between the separatrices: in the asymp- 
totic limit $ — > 0, u)(E) approaches a rectangular form. 
Besides, the amplitude q n of the nth Fourier harmonic 
of the oscillation asymptotically approaches a non-small 
constant value in the whole energy range between sepa- 
ratrices. These two properties are responsible for most 
of the characteristic features of the global chaos onset 
in between the separatrices. The most striking one is a 
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drastic facilitation of the global chaos onset when the per- 
turbation frequency ujf approaches uj m or its multiples: 
the perturbation amplitude h gc required for global chaos 
possesses, as a function of the perturbation frequency uif, 
deep spikes close to uj m or its multiples. 

On the basis of the theory for the boundaries of the 
chaotic layers, we have developed a self-consistent asymp- 
totic theory for the spikes in the vicinity of the min- 
ima. In particular, the explicit asymptotic expressions 
for the very minima are given in Eqs. (52) and (55). 
The minimal amplitude h gc is smaller than the typical 
h gc for ujf beyond the close vicinity of u> m by a factor 
~ 10R ~ 101n(l/$). The theory well agrees with the 
simulation results. 

We have also found the mechanisms responsible for 
the spike wings (Figs. 9, 11). The theory well fits the 
simulations (Fig. 10). The asymptotic shape of the spike 
is universal: it is described by Eq. (67) (Fig. 12). 

The facilitation of the global chaos onset may have the 
following applications in particular: 

• drastic increase of the DC-conductivity of a 2D elec- 
tron gas in a ID magnetic superlattice; 

• significant decrease of the activation barrier for 
noise-induced escape over double/multi-barrier 
structures, that may lead to a drastic acceleration 
of the diffusion in periodic structures; 

• strong facilitation of the stochastic web formation. 
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APPENDIX A: SEPARATRIX MAP ANALYSIS 

The chaotic layers of the system (2) associated with the 
separatrices of the unperturbed system (1) are described 
here by means of the separatrix map. To derive the map, 
we follow the method described in but the analysis 
of the map significantly differs from existing ones 0, 0, 
EL H ■ It constitutes the method allowing to calculate 
chaotic layer boundaries in the phase space (rather than 
just in energy). It also allows to quantitatively describe 
the transport within the layer in a manner different from 
existing ones (cf. [H, [6(| and references therein) . 

We present below a detailed consideration of the lower 
chaotic layer while the upper layer may be considered 
similarly and we present just the results for it. 




time, 



FIG. 14: Schematic example of the time dependence of the 
velocity of the perturbed system (thick solid line) in the case 
when the energy of motion varies in the close vicinity of the 
top of the lower potential barrier. The dashed line marks the 
zero level of the velocity. Pulses of the velocity are schemati- 
cally singled out by the parallelograms (drawn by a thin solid 
line). The two sequences of time instants (..., ti-i, U, U+i, — ) 
and (...,t' i _ 1 ,t' i ,t' i+1 , ...) correspond to beginnings and centers 
of the pulses, respectively. 



1. Lower chaotic layer 



a. Separatrix map 

The typical function q(t) for the trajectory close to 
the inner separatrix (the separatrix corresponding to the 
lower potential barrier) is shown in Fig. 14. One can 
resolve pulses in q(t). Each of them consists of two ap- 



proximately antisymmetric spikes [611 ]. The pulses are 
separated by intervals during which \q\ is relatively small. 
Generally speaking, successive intervals differ of each 
other. Let us introduce the pair of variables E and ip: 



H , ip = Uft + ip a . 



(Al) 



where the constant Lp a may be chosen arbitrarily. 

The energy E changes only during the pulses of q(t) 
and remains nearly unchanged during the intervals be- 
tween the pulses, when \q(t)\ is small Q. We assign 
numbers i to the pulses and introduce the sequences of 
(Ei,(fi) corresponding to initial instants of pulses U. In 
such a way, we obtain the following map (cf. Q): 



E i+1 =E i +AE i , (A2) 
w /7 r(3 - sign(£ i+1 - £ 6 (1) )) 



tfii+i = <Pi + 
AE, = h 



2u(E, 



dt q(t) cos(u)ft), 
ith pulse 



where J lt j 1 pu i se means integration over the ith pulse. 
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Before deriving a more explicit expression for AEi , we 
make two following remarks. 

1. Let us denote with t\ the instant within the ith 
pulse when q is equal to zero (Fig. 14). The function 
q(t — i'j) is an odd function within the ith pulse and it 
is convenient to transform the cosine in the integrand in 
AE % (A2) as 



cos(ujft) = cos(w/(i — t'j) + Uft'i) = 

cos(w/(f - t'i)) 008(^/4) - sin(cj/(i - *•)) sin(w/0> 

and to put ip a = — ti), so that ipi = Wft^. 

2. Each pulse of q contains one positive and one nega- 
tive spike. The first spike can be either positive or nega- 
tive. If E changes during the given nth pulse so that its 
value at the end of the pulse is smaller than E^ , then 
the first spikes of the ith and (i + l)st pulses have the 
same signs. On the contrary, if E at the end of the ith 
pulse is larger than E^ , then the first spikes of the ith 
and (i + l)st pulses have opposite signs. Note that Fig. 14 
corresponds to the case when the energy remains above 
E^ during the whole interval shown in the figure. This 
obviously affects the sign of AEi , and it may be explicitly 
accounted for in the map if we introduce a new discrete 
variable (jj = ±1 which characterizes the sign of q at the 
beginning of a given ith pulse, 



(Ti = sign(g(£i)) , 
and changes from pulse to pulse as 



(A3) 



= en sign(^ 1} - E i+1 ) . (A4) 

With account taken of the above remarks, we can rewrite 
the map (A2) as follows: 



E l+1 = Ei + vihe^smiwi 



<Pi+i 

0~i+l 

Mow) 



'-Pi 



q; / 7r(3-sign(£ i+1 -£f ) )) 



Ui sign(^ 1} - E i+1 ), 



Mow) 



- <;,; / dt q(t - t'j) sm(u f(t - i-)) 

ith pulse 

rU+i 

-loi J dtq(t-t$sw(uf(t-ti)). 



The map similar to (A5) was introduced for the first 
time in [62| , and it is often called as the Zaslavsky sepa- 
ratrix map. Its mathematically rigorous derivation may 
be found e.g. in the recent major mathematical review 
[32| . The latter review describes also generalizations of 
the Zaslavsky map as well as other types of separatrix 



maps. The analysis presented below relates immediately 
to the Zaslavsky map but it is hoped to be possible to 
generalize it for other types of the separatrix maps too. 

The variable e( /otu ) introduced in (A5) will be conve- 
nient for the further calculations since it does not depend 
on i in the lowest-order approximation. A quantity like 
Si = ft,|e(' ou, )| is sometimes called the separatrix split Q 
since it is conventionally assumed that the maximal de- 
viation of energy on the chaotic trajectory from the sep- 
aratrix energy is of the order of Si @, H, 0, [f| . Though we 
shall also use this term, we emphasize that the maximal 
deviation may be much larger. 

Dynamical chaos appears in the separatrix map (A5) 
because uj(E — > E^ 1 ') — > 0. Various heuristic criteria 
were suggested for the estimate of the chaotic layer width 
in energy 0, H, H, @ . Frequencies relevant to our problem 
are much smaller than the reciprocal width of the spikes 
of q(t). For such frequencies, all these criteria [1, [3,13, 
give: 



UJQ 



.(low) 



Uf < U> W 1, 



(A6) 



where ajg is the frequency of eigenoscillation at the bot- 
tom of the potential well. 

The estimate (A6) was used in our earlier theory [l2l |. 
But we found later that, for the case of small ujf, the 
aforementioned criteria were insufficient, so that the es- 
timate (A6) was incorrect [63[ (cf. also [13, [fill). More- 
over, to search a uniform width of the layer is incorrect 
in cases like ours, where the width strongly depends on 
the angle. At the same time, the lowest- order formulas 
for the spike minimum (hs,co s ) are not affected by this, 
so that the results of 12] (with only the lowest-order 
formulas) are correct. Still, the higher-order corrections 
(quite significant for hi^ if $ is moderately small) would 
be incorrect if they were calculated on the basis of the 
estimate (A6). Besides, the paper [l2[ did not address 
the intriguing question: why does even a small excess 
(A5) OI h over h gc (ujf) result in the onset of chaos in a large 
part of the phase space between the separatrices, despite 
the fact that the width of the chaotic layers associated 
with the nonlinear resonances is exponentially small for 
h = h gc (ujf)? The analysis of the separatrix map pre- 
sented below resolves these important problems. 

In the adiabatic limit u>f — > 0, the excess of the upper 
boundary of the lower layer over the lower barrier 

also 



Ef does not depend on angle and equals 2irh [63| (cf. 

[Hj|). But tof relevant for the spike of h gc (ujf) can- 
not be considered as an adiabatic frequency, despite its 
smallness, because it is close to uo m or to its multiple 
while all energies at the boundary lie in the range where 
the eigenfrequency is also close to Lu m : 



3 = 



*(2i- 
1,2,3, 



l)u) m « (2j 



i)^(4 1} ) 



(A7) 
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The validity of (A7) (confirmed by the results) is cru- 
cial for the description of the layer boundary in the rel- 
evant case. 

b. Separatrix split 

Let us explicitly evaluate ^ low \ Given that the energy 
is close to the velocity g(t-^) in t {low) (A5) may be 
replaced by the corresponding velocity along the separa- 
trix associated with the lower barrier, qs° w \t — while 
the upper limit in the integral may be replaced by infin- 
ity. Besides, in the asymptotic limit $ — > 0, the interval 
between spikes within the pulse becomes infinitely long 
[6l| and, therefore, only short (~ i-^ 1 ) intervals corre- 
sponding to the spikes contribute to the integral in g( /ou ') 
(A5). In the scale wj 1 , they may be considered just as 
two "instants": 

fg/ 2) -^±-^, $-0. (A8) 

In the definition of e'- low ^ (A5), we substitute the argu- 
ment of the sine by the corresponding constants for the 
positive and negative spikes respectively: 

\ 4W m / Jp 0s itivc spike 

<A9) 

$ -> 0. 

In the derivation of the first equality in (A9), we have also 
taken into account that the function q£ ow ' (x) is odd. In 
the derivation of the second equality in (A9), we have 
taken into account that the right turning point of the 
relevant separatrix is the top of the lower barrier and the 
distance between this point and the left turning point of 
the separatrix approaches ir in the limit $ — > 0. 

For the frequencies relevant to the minima of the spikes 
of h gc (ujf), i.e. for u>f = u>s rs (2j — l)w m , we obtain: 

£ (*<™)( w 0)) ra 2tt sin ((2j - 1) J) = VM-l)^], 
j = 1,2,3,..., $^0. (A10) 

For moderately small $, it is better to use the more 
accurate formula: 

/>CO 

e {low) {cj f ) = 2 / dt qV° w) (t)sm{LU f t), (All) 
Jo 

where the instant t — corresponds to the turning point 
of the separatrix to the left from the lower barrier, i.e. 
q { J ow \t = 0) = while qi° w) > for all t > 0. The 
dependence |e^ 0, "'(w/)| by Eq. (All) is shown for $ = 
0.2 in Fig. 15(a). For small frequencies, the asymptotic 
formula (A9) well fits the formula (All). 




FIG. 15: The theoretical estimates for the normalized sep- 
aratrix split (for $ = 0.2) as a function of the perturbation 
frequency, for the lower and upper layers in (a) and (b) re- 
spectively. The solid lines are calculated by Eqs. (All) and 
(A43) (for (a) and (b) respectively) while dashed lines repre- 
sent asymptotic expressions (A9) and (A44) respectively. 



c. Dynamics of the map 

Consider the dynamics of the map (A5), when uif is 
close to the spikes minima: u>f « nu> m where n = 2j — 1 
while j — 1, 2, 3, . . .. Let the energy at the step i = — 1 be 
equal to . The trajectory passing through the state 
with this energy is chaotic since (w(-E)) -1 diverges as 
E — > Ejp and, therefore, the angle ip-i is not correlated 
with the angle on the previous step (p-2- The quantity 
er_i is not correlated with er„2 either. Thus, sin(</?_i) 
may take any value in the range [—1, 1] and cr_i may 
equally take the values 1 or -1. Therefore, the energy 
may change on the next step by an arbitrary value in the 
interval [-h\e^ low \ h\e ( - low '>\]. Thus, E - E^ may have 
a positive value ~ h\e^ low ^\ Then, the approximate 
equality juo(Eq) ~ u> m holds, provided the value of h is 
from the relevant range. Allowing for this and recalling 
that we are interested only in those realizations of the 
map such that Eq > E^\ the relevant realization of the 
map i = — 1 — * i = may be written as: 

E = E^ +a_i/ie^sin(^_ 1 ) = 

= E^+ftle^sin^i)!, 
ipo w ip-! + mr, 

<T =-<7_l. (A12) 
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One may expect that the further evolution of the map 
will approximately follow, for some time, the trajectory 
of the system (4) with the initial energy Eo (A12), and 
with an arbitrary p-\ and initial slow angle ip somehow 
related to tpo ~ P-i + nix . Let us prove this explicitly. 

Consider two subsequent iterations of the map (A5): 
2i — > 2i + 1 and 2% + 1 — * 2i + 2 with an arbitrary 
% > 0. While doing this, we shall assume the valid- 
ity of (A7) (it will be clarified below when this is true) 
from which it follows that: (i) u(Ek+i) ~ uj(Ek), (ii) 
Pk+i — <pk ~ nn = (2j — l)ir. It will follow from the 
results that the neglected corrections are small in com- 
parison with the characteristic scales of the variation of 
E and ip (cf. the conventional treatment of the nonlin- 
ear resonance dynamics 0, S S, 0, 0, 0). Besides ' [t 
follows from (A5) that, while the energy remains above 
the barrier energy, cj, oscillates, so that a 2i — Co an d 
C2i+i = -co- Then, 



E21+1 

P2i+1 



E 2l + aohe^ ow hm( V2t ), 



^2 



w{E 2i +i) 

UJf - naj(E 2 i) 



<P2i + nir + n 



lu(E : 



2i J 



(A13) 



E- 



2i+2 



P2i+2 



a he^ low> sm((p 2t+ i) 



Low) j 
Low) 



= E 2 i+i 

= E 2l+1 + a he {low} sin(^ 2l +i - nir) 
w E 2 i + ao2h^ low hm(ip2i), 



f2i+l 



(f 2 i + 2ml + 2ir 



uj(E 2i+2 ) 

UJf - nu)(E 2i ) 



2* I 



(A14) 



(the second equality in the map for E 2 i +2 takes into ac- 
count that n is odd so that sin((^ — mi) = — sin(</?)). 

The quantity p2%+2 — pn — 27m is small, so the map 
2i — > 2i + 2 (A14) may be approximated by differential 
equations for E 2 i and tp 2 i = p 2 i — 2imi: 



dE 2l 
d(2i) 

dlfi 2 j _ 7T 

d(2i) ~ uj(E 2i ) 
P2i = P2i - 2nni 



he«° w Um(<p2i), 



(uj f - nuj(E 2i )), 



(A15) 



Let us (i) use for e^ lowS > the asymptotic formula (A10), 
(ii) take into account that the increase of i by 1 corre- 
sponds to the increase of time by n/u>(E), and (iii) trans- 
form from the variables (E, p) to the variables (/, ip = 
mr(l — ao)/2 — (p). Equations (A15) reduce then to: 



dl 

dT 



-hV2(-l)i^] sin(^), 



<hjj_ 



nu! — ujf, 

1 - (T 



(A16) 



n = 2j - 1. 



Equations (A16) are identical to the equations of mo- 
tion of the system (4) in the lowest-order approximation, 
i.e. to the equations (14) where q n is replaced by its 
asymptotic value (12) and the last term in the right-hand 
part of the second equation is neglected, being of higher 
order in comparison with the term nui — Wf. 

Apart from the formal identity of Eqs. (A16) and (14), 
tp in (A16) and ip in (14) are identical to each other. 
Necessarily t\ corresponds to a turning point (see Fig. 
14) while the corresponding ip is equal to 2ni or n + 2i:i 
for the right and left turning point respectively (see (4)) 
i.e. ip = 2ni + tt(1 — <7j)/2, so that '0(14) = wl 3 ~ u) f't — 
tmt(1 - a)/ 2 -Cp = ij)(Axg). 

The relevant initial conditions for (A16) follow from 
(A12) and from the relation between ip and p: 



1(0) = I(E = E^ + hV2n\ sin(^(0))|), (A17) 

while ip(0) = ?i7r(l — (7o)/2 — </3o m ay be an arbitrary angle 
from the ranges where 



-1)["/ 4 1 sin(^(0)) < 0. 



(A18) 



For moderately small <&, it is better to use the more 
accurate dynamic equations (14) instead of (A16) and the 
more accurate initial value of action instead of (A17): 



1(0) = I(E = E, 



(i) 



01 1 sin 



mom, 



(A19) 

with e {low ^ calculated by (All). 

We name the quantity <5/| sin('0)| the generalized sepa- 
ratrix split (GSS) for the lower layer. Unlike the conven- 
tional separatrix split Si 4], it is angle- dependent. The 

curve I(rp) = I(E = E^ + Si\sin(ip)\) may be called 
then the GSS curve for the lower barrier and denoted as 

Finally, let us investigate an important issue: whether 
the transformation from the discrete separatrix map (i.e. 
(A13) and (A14)) to the differential equations (A15) is 
valid for the very first step and, if it is so, for how long 
it is valid after that. The transformation is valid as long 
as u>(Ek) ~ nu>f i.e. as long as E^ is not too close to 

the barrier energy E^\ At the step k — 0, the system 
stays at the GSS curve, with a given (random) angle ^(0) 
from the range (A18). Thus, at this stage, the relation 
(A7) is certainly valid (for the relevant range of h and 
for an angle a little away from the narrow vicinity of the 
multiples of 7r) . 

The change of energy at the next step is positive too: 
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Et-Eo = aohe^ sin(^ ) « 

w -a^he {low) sin(^_i - utt) = 

= o-_i/ie^ sin(^-i) =E - £7_i > 0. 

This may also be interpreted as a consequence of the first 
equation in (A16) and of the inequality (A18). 

Hence, (A7) is valid at the step k = 1 too. Similarly, 
one can show that E% — E\ > 0, etc. Thus, the trans- 
formation (A13,A14)— >(A15) is valid at this initial stage 
indeed, and the evolution of (E, (p) does reduce to the 
resonant trajectory (14) with an initial angle from the 
range (A18) and the initial action (A19). This lasts until 
the resonant trajectory meets the GSS curve in the adja- 
cent 7r range of -0 i.e. at t such that the state (I(t), ip(t)) 
satisfies the conditions: 

m = 4ss(#))> imM - 0(0)/*] = I- (A20) 

At this instant, the absolute value of the change of energy 
Ek in the separatrix map (A13) is equal to Ek — E^ 1 ' 
(just because the state belongs to the GSS curve) but the 
sign of this change is negative since the sign of sin(tpk) 
is opposite to the sign of sin(</?o)- Therefore, at the step 
k + 1, the system gets to the very separatrix, and the 
regular-like evolution stops: at the next step of the map, 
the system may either again get to the GSS curve with 
a new (random) angle from the range (A18) and start a 
new regular-like evolution as described above; or it may 
get to the similar GSS curve below the barrier and start 
an analogous regular-like evolution in the energy range 
below the barrier, until it stops in the same manner as 
described above, etc. 

This approach makes it possible to describe all features 
of the transport within the chaotic layer. Their detailed 
description will be done elsewhere while, in the present 
context, it is most important to describe the upper outer 
boundary of the layer. 



d. Boundary of the layer 

We may now analyze the evolution of the boundary of 
the layer as h grows. Some of the evolution stages are 
illustrated by Figs. 8, 9 and 16. 

It follows from the analysis carried out in the previous 
subsection that any state (in the I — ip plane) lying be- 
yond the GSS curve but belonging to any trajectory fol- 
lowing the equations (14) which possesses common points 
with the GSS curve belongs to the chaotic layer: the sys- 
tem starting from such a state will get, sooner or later, 
to the separatrix, where the chaotization will necessar- 
ily occur. Therefore, the upper boundary of the chaotic 
layer coincides with the trajectory following the equa- 
tions (14) with the initial action (A19) and an initial 




angle ij/ 

FIG. 16: A schematic figure illustrating the formation of the 
boundary of the lower chaotic layer for h < hcr(tJf) in the 
ranges of relevant to (a) odd, and (b) even spikes. The 
dashed magenta line shows the GSS curve in the energy-angle 
plane: E{^) = E { $ ss $) = E^ + 5i\sui{$)\. Green lines 
show examples of those trajectories (14) which have points in 
common with the GSS curve. One of such trajectories (14) 
(shown by the thick green line) forms the upper boundary of 
the lower chaotic layer: in (a), it is the trajectory tanyent to 
the GSS curve; in (b), it is the upper part of the separatrix 
generated by the saddle "s" . Yellow dots indicate the relevant 
common points of the GSS curve and the thick green line. 
They have angles ±V>t and energy E\ in the case (a), and 
angles and energy E^ in the case (b). The minimum 

and maximum deviation of energy on the boundary from the 
barrier energy are denoted as S^ ln and 5max respectively. The 
maximum deviation on the GSS curve is equal to Si. 



angle ip(0) from the range (A18) such that the trajec- 
tory deviates from the barrier energy more than a trajec- 
tory (14)-(A18)-(A19) with any other initial angle does. 
There may be only two topologically different options for 
such a trajectory: either it is the trajectory tangent to 
the GSS curve, or it is the separatrix trajectory which 
intersects the GSS curve (some schematic examples are 
shown in Figs. 16(a) and 16(b) respectively; some real 
calculations are shown in Figs. 8 and 9). 

1. Relatively small h 

Consider first values of h which are large enough for the 
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condition (A7) to be satisfied (the explicit criterion will 
be given in (A31)) but which are smaller than the value 



he) = hcr(uf) determined by Eq. (41) (its meaning is 
explained below). The further analysis within this range 
of h differs for the ranges of uif relevant to odd and even 
spikes, and we consider them separately. 

A. Odd spikes 

The relevant frequencies are: 



-,(') 



uj f w noj m , 
3 = 1)3,5, . 



2j - 1, 



(A21) 



Let us seek the state {I^j-tpP} (with within the 
range ]0, 7r[) where the resonant trajectory curve is tan- 
gent to the GSS curve. With this aim, we equal both the 
actions and the derivatives of both curves. The equality 

4 l) = HE = 

The derivative along the GSS curve 



of actions immediately yields 1^ via 



E 



(0 



7, 



GSS 



is obtained by differentiation of 7q SS (?/>). The derivative 
along a resonant trajectory can be found dividing the 
first dynamic equation in (14) by the second one. Sub- 
stituting the expression of if 1 via ipt into the equality 
of the derivatives, we obtain a closed equation for ipP , 
and its solution immediately gives us the relevant 4>(0): 



E (,o,B) |cas($ ) 
+q n (E) sin($°) 



E=E 



m 



nu(E) 
0, 



dE 



(A22) 



E 



E, 



(i) 



+ / l |e^)|sin(#), 
n = 2j - 1, j 



& W G [0,7T], 



1,3,5,. 



4,(0) 



7(i) 

= 4>i 



A careful analysis of the phase space structure shows 
that, in the present case (i.e. when h < hf£}(u>f) while j 
is odd), there is no separatrix of the resonant Hamilto- 
nian (4) which would both intersect the GSS curve and 
possess points above the tangent trajectory (6?]]. Thus, 
for this range of h, the outer boundary of the chaotic 
layer is formed by the trajectory following the dynamical 
equations (14) with the initial angle by (A22) and initial 
action (A19) (Fig. 16(a)). 

Let us find the lowest-order solution of Eq. (A22). We 
neglect the term 1 — ujf/(nuj(E)) (the result will confirm 
the validity of this) and use the lowest-order expression 
for the relevant quantities: namely, Eqs. (A10) and (12) 
for {( low ) and q n respectively, and the lowest-order ex- 
pression for dq n /dE which can be derived from Eq. (11): 



dgn(g) 

dE 
n = 2j- 1, 



[-1)1* 



4^2 [E - 75 b (1) ) In (S" 1 ) 



E - El L> < $ -> 0. 



(A23) 



Then Eq. (A22) reduces to the following equation 



tan (Vt ) - 81n 



(A24) 



The lowest-order solution of (A24) in the range ]0,7r[ 
reads as 



$° = ("l) [f] 



81n(l/$) 



i - (-i)W 



(A25) 



It follows from the definition e[ 1) (A22) and from (A25) 
that the lowest-order expression for Ejp — 25^ reads as 



r 3/2 



\fiMT/W) In 



(A26) 



The next step is to find the minimal value of the energy 
on the boundary of the layer, E^ n . It follows from the 
analysis of the dynamical equations (14) that the cor- 
responding angle V'min is equal to if sign(q2j-i) > 
(i.e. j = 1,5,9, . . .) or to tt if sign(g 2j _i) < (i.e. j = 
3,7,11,...): cf. Fig. 8(a). Given that the Hamiltonian 
(4) is constant along any trajectory (14) while the bound- 
ary coincides with one of such trajectories, the values of 
the Hamiltonian (4) in the states {I(E^ in ), tp = ip m in} 
and {It } should be equal to each other. In the 
explicit form, this equality may be written as 



(0 



dE 1- 



nw(E) 



(A27) 



Let us find the lowest-order solution of Eq. (A27). As- 
sume that -E™L still belongs to the range of E where 
w(E) « uj m (the result will confirm this assumption). 
Then the integrand in (A27) goes to zero in the asymp- 
totic limit $ — * and, hence, the integral may be 
neglected (the result will confirm the validity of this). 
The remaining terms in Eq. (A27) should be treated 
very carefully. In particular, it is insufficient to use 
the lowest-order value (12) for q n since it is the differ- 



ence between q n (E*jp) and q n {E^^) that matters. More- 
over, the approximate equality q n (EP) — qn(E^ in ) 



dqM l) )/^Ef\E[ l) 



T^min) does not apply here ei- 



ther since, as it follows from Eq. (A23), the derivative 
dq n (E)/dE may strongly vary in the range [E^ in , E^] 



(0 



E^ ') > 1 (the result will show 



that it is the case) . That is why it is necessary to use for 
q n the more accurate expression (11). Allowing for the 
asymptotic expression (A25) of ?/>j and keeping only the 
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lowest-order terms, one can finally reduce Eq. (A27) to 
the relation 



At the same time, unlike in the case of odd spikes, 
there exists a saddle with an angle 



In 



E 



E, 



(i) 



, ^min _ E b 



(1) 



(A28) 



Substituting here the asymptotic value of (A26), we 
obtain the final lowest-order expression for the minimal 
(along the boundary) deviation of the energy from the 
barrier: 



(0 _ 



.l-(-l)L 



(A33) 



while the energy (which may be found as the appropriate 
solution of Eq. (15)) lies in the relevant vicinity of the 
lower barrier (Fig. 16(b)). In the lowest-order approxi- 
mation, this saddle energy reads: 



min 



^3/2 



E y ' - E 



(1) 



(El 



(A29) 



2eVa y/]n(l/$) /n 

It is necessary and sufficient that the condition u>(E) w 
Lj m is satisfied at the minimal and maximal energies of 
the boundary to ensure that the second equality in (A7) 
holds true, i.e. that u>(E) is close to u> m for all points of 
the boundary. 

At the minimal energy, this condition reads as 



<5min) « W» 



(A30) 



Eq. (A30) determines the lower limit of the relevant 
range of h. The asymptotic form of (A30) is: 



2V2ln(ln(4e/$)) 



(A34) 



This saddle (denoted in Fig. 16(b) as V) generates a 
separatrix. Its upper whiskers go to the similar adjacent 
saddles (shifted in ip by 2n). In the asymptotic limit 
$ — > 0, the upper whiskers are much steeper than the 
GSS curve and hence they do not intersect it [6j|. As 
concerns the lower whiskers, they do intersect the GSS 
curve and, moreover, two intersections lie in the relevant 
energy range (Fig. 16(b)). Let us show this explicitly. 
Let us write the expression for the Hamiltonian (4) in the 
relevant vicinity of the barrier energy (i.e. where Lu m — 
w{E) <§; w m ), keeping, in the expression, both the lowest- 
order terms and the terms of next order (in particular, 
we use Eq. (11) for q n (E) and take into account that 
< V2 — nq n (E) <C \/2 for the relevant range of E) : 



In 



^Y/inCi/*) 



ln(l/$) 



< 1. 



(A31) 



We emphasize that any h of the order of h s o (52) satisfies 
this condition. In the asymptotic limit <E> — > 0, the left- 
hand part of Eq. (A31) goes to zero. 

As for the maximal energy, it may take values up to the 
energy of the lower saddle "si", i.e. E s i (18). Obviously, 
(A7) is valid at this saddle, too. 

B. Even spikes 

The relevant frequencies are: 



ujf w nuj m , 
J = 2, 4, 6,. 



n = 2j - 1, 



(A32) 



In this case, q n {E) and dq n (E)/dE have different signs 
for all E within the relevant range (i.e. where uj(E) w 
Wm, q n (E) w q n (E m )): cf. (12) and (A23). Then, in 
the asymptotic limit $ —> 0, Eq. (A22) for the tangency 

does not have any solution for ip^ in the relevant range 
[68l |. There may be only solutions very close to some 

of 7r integers, and the corresponding energies E^ are 
very close to E^ i.e. ui(E^) <C ui m : therefore they are 
irrelevant. 



H(I = I(E = E[ 



21n(f) 



uj m — uj(E + S) <C uj„ 




(A35) 



The Hamiltonian H should possess equal values at the 
saddle "s" and at the intersections of the separatrix and 
the GSS curve. Let us denote the angle of the intersection 
in the range ]0, 7r[ as ipf', and let us denote the deviation 
of its energy E^ from E^ as sf^ = Si sin(^P). 

Assuming that I'tpf^ — ips^\ *C 1 (the result will confirm 
this) so that cos^f) » (-l)["/ 4 l(l-(^f -^°) 2 /2) w 
(_l)[«/4](i _ $ l) /5i) 2 /2) « J-l)[-/4] (1 _ (5(0/^)2/4), 
the equality of the values of H reads as: 



21n(f) 



= hV2 



In [4h 



(S 



(0\2 



2^2h ' 



(A36) 
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Let us assume that, in the asymptotic limit $ — > 0, 
<C <5 S (the result will confirm this). Then the left- 
hand part is asymptotically smaller than the first term 
in the right-hand part. So, Eq. (A36) implies, in the 
asymptotic limit, that the right-hand side equals zero. 
Expressing h via 6 S from Eq. (A34), we finally obtain a 

closed transcendental equation for 5 a /8f : 




n (In (In (#)))< 



= A. 



(A37) 



In the asymptotic limit $ — > 0, the quantity A diverges 
and, hence, the lowest-order asymptotic solution of Eq. 
(A37) reads as 



5s_ 



I 2A 
\MA)' 



(A38) 



Substituting here the expression (A34) for 5 S and the 
expression (A37) for A, we obtain: 



s(0 



ln(ln(f)) 



In 



(A39) 



Thus, we have proved the following asymptotic prop- 
erties of the separatrix generated by the saddle "s": 1) 
the lower whiskers of the separatrix do intersect the GSS 
curve in the relevant range of E (i.e. where the res- 
onant approximation is valid), 2) the upper whiskers of 
the separatrix do not intersect the GSS curve (there is no 
solution of Eq. (A36) in the range tfp > S s ). The former 
property confirms the self-consistence of the asymptotic 
theory for even spikes; the latter property means that the 
upper outer boundary of the lower chaotic layer is formed 
by the upper whiskers of the separatrix generated by the 
saddle "s". 

Finally, we explicitly note that the minimal (along the 
boundary) deviation of energy from the barrier energy 
occurs occurs right at the saddle "s" , i.e. 



trajectory starting from the state with the angle (A22) 
and action (A19), for odd spikes, or from the saddle "s", 
for even spikes, is encompassed by the resonance separa- 
trix. So, it does not form the outer boundary of the layer 
anymore. Rather it forms the inner boundary i.e. the 
boundary of the main island of the stability inside the 
layer, repeated periodically in tp with a period 2ir (cf. 
analogous islands in the upper layer in Fig. 8). Unless 
the lower chaotic layer reconnects with the upper one, the 
outer boundary of the lower layer is formed by the upper 
part of the resonance separatrix. The relevant initial an- 
gle i/)(0) on the GSS curve corresponds to the intersection 
of the GSS curve with the resonance separatrix (cf. the 
analogous situation for the upper layer in Fig. 8). 



2. Upper chaotic layer 

The upper chaotic layer may be treated analogously 
[70l ] to the lower layer. We present here only the results. 

Similarly to the lower-layer case, one may consider 
the ranges of relatively small h (namely, smaller than 
/icr = h^\u)f) determined by Eq. (42)) and relatively 

large h (i.e. h > h^fi). In the former range, the forma- 
tion of the boundary occurs in a manner which is, in a 
sense, opposite to that for the lower-layer case. For even 
spikes, the lower outer boundary is formed by tangency 
while, for odd spikes, it is formed by the lower part of 
the separatrix generated by the saddle "s" , analogous to 
the saddle "s"in the lower-layer case pZll j - 



So, for even spikes, the angle of tangency ipt 
termined by the following equation 



(«) 



is de- 



i e ^icos(^r ; )(i 



-q n (E) S in(ti u) ) 



E=El 



nuj(E) 
= 0, 



dE 



(A41) 



E 



(«) = E (2) 



^ (u) e[0,7r], 



and ifil determines the tangency energy: 



, -h\e^\ S m(4 u) ) 
n = 2j- 1, 

(«) 



.7 = 2,4,6, 



5^ = 8- (A40) 

2. Relatively large h. 

As h grows, the boundary of the layer raises up while 
the lower part of the resonance separatrix, on the con- 
trary, goes down. They reconnect at the critical value 
of h, her = hcr(u>f), determined by Eq. (41), that may 
be considered as the absorption of the resonance by the 
chaotic layer. If h grows further, then the GSS curve 
and the resonance separatrix intersect. As a result, the 



eW=4 2 > -h\e^\sm($ u) ), (A42) 
where the quantity e^ up ^ is described by the formula 

e (" p ) {u>f)=2 dt (t) cos(w ft) , (A43) 

where qi up \t) is the time dependence of the velocity 
along the separatrix associated with the upper barrier 
and the instant t = is chosen so that qi up \t = 0) 



28 



is equal to the coordinate of the lower barrier while 
qi up) > for t £ [0,oo[. The dependence |e<»(u;/)| 
in Eq. (A43) is shown for $ = 0.2 in Fig. 15(b). 
The asymptotic form of Eq. (A43) reads as 

e (up) _ e («p)( ) = 27rcos [ ™/ ] _ ( A44 ) 

For w/ = wr « (2j - l)w m , Eq. (A43) reduces to 

e (up) (t^) « 27rcos((2j - = \/27r(-l)[^], 

.7 = 1,2,3,..., $^0. (A45) 

The lowest-order solution of (A41) is given in Eq. (37), 

so that is approximated by Eq. (38). The maximal 
energy on the lower boundary of the layer corresponds 



to $(t) = 7r if j = 2, 6, 10, ... or if j = 4, 8, 12, . . . and 
is determined by Eq. (39). The asymptotic value of the 
minimal deviation from the upper barrier of the energy 
at the boundary, 5^ n , is given in Eq. (40). 

For odd spikes, the boundary is formed by the lower 
part of the separatrix generated by the saddle "s " . The 
angle of the saddle is given in Eq. (33) while the deviation 
of its energy from the barrier is approximated, to the 
lowest-order approximation, by Eq. (34). 

As h grows, the boundary of the layer goes down while 
the upper part of the upper resonance separatrix goes up. 
They reconnect at h = n£/ = /i^ (w/), as determined 
by Eq. (42) , that may be considered as the absorption of 
the resonance by the chaotic layer. 

For larger h, the boundary of the layer is formed by 
the lower part of the upper resonance separatrix (Fig. 8), 
unless the latter intersects the lower GSS curve (in the 

latter CclS6, Her marks the global chaos onset). 
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